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SPACECRAFT  ATTITUDE  CONTROL 
FOR  EXTENDED  MISSIONS 


ABSTRACT 

This  thesis  considers  the  problem  of  providing  attitude 
control  for  a  spacecraft  engaged  in  an  extended  mission.     As  a 
basis  for  the  choice  of  a  suitable  attitude  control  system  the 
following  requirements  are  applied. 

Maximum  reliability 
Minimum  ejection  of  mass 
Minimum  average  power 
Minimum  system  weight 
Minimum  peak  power 

An  interplanetary  mission  of  400  days  duration  is  adopt- 
ed as  a  general  guide  for  the  problem,  but  most  of  the  equations 
and  comparisons  are  presented  in  parametric  form.     Extended 
missions  imply  that  a  momentum  exchange  type  attitude  control 
system  be  used  to  minimize  ejection  of  fuel  mass,   and  the  thesis 
primarily  considers  only  systems  of  this  type.     The  thesis  derives 
the  equations  of  motion  for  a  spacecraft  equipped  with  eighteen 
different  control  systems.     The  control  system  chosen  to  best 
satisfy  the  five  design  requirements  is  a  system  consisting  of 
four  gyro-type  controllers  arranged  in  two  pairs  with  each  pair 
operating  back-to-back  to  minimize  control  cross  coupling 
torques.     One  pair  of  controllers  provide  roll  torques,   the  other 
pair  provides  pitch  torques,   and  all  four  controllers  contribute 
yaw  torques. 

The  stability  and  control  analysis  considers  operation 
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contribute  yaw  torques. 

The  stability  and  control  analysis  considers  operation 
of  the  spacecraft  in  three  modes. 

Zero  Input  Mode 
Rate  Control  Mode 
Position  Control  Mode 

Each  of  the  modes  are  evaluated  for  roll  motion  by  assunning 
negligible  interaxial  coupling,   and  the  analysis  includes  opera- 
tion of  the  controller  gimbal  angles  to  large  angles. 
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CHAPTER  6 
CLOSING  THE  ATTITUDE  CONTROL  LOOP 

6.  1 Introduction 

Previous  chapters  have  defined  the  physical  character- 
istics of  the  spacecraft,   have  chosen  the  four  gyro  (12-34-1234) 
system  as  being  the  most  attractive  control  system,   and  have 
defined  probable  disturbances  on  the  spacecraft.     Before  the 
mathematical  models  representing  the  spacecraft  and  control 
system  can  be  combined  into  a  single  loop  it  is  necessary  to 
specify  an  attitude  sensor  and  some  type  of  torque  motor  which 
will  receive  signal  information  from  the  sensor  and  provide  the 
proper  angular  rates  and  positions  to  tlie  spin  axes  of  the  four 
controllers.     The  sequence  of  this  chapter  follows  an  order  that 
first  determines  the  dynamics  of  the  individual  controllers  by 
methods  similar  to  those  of  Chapter  2  which  considered  the  com- 
plete spacecraft.     Next,   the  motors  and  drive  systems  which 
drive  the  controllers  are  defined.     Then,   the  sensor  and  related 
amplifiers  are  assumed  together  with  the  provision  for  tandem 
compensation,   if  required.     Finally,   all  of  the  above  components 
are  combined  into  a  single  loop. 

6.  2  Response  of  the  Controllers  to  Moments  Applied  by 

the  Roll  and  Pitch  Torque  Generators 


In  a  manner  parallel  to  that  of  Section  2.  3  the  funda- 
mental equation  of  motion  of  a  controller  consisting  of  a  gyro 
case,   a  gimbal,   and  a  rotor  is 

-4H-      =  y  M  ,  (Eq  6.2.  1) 

dtj  J       £_,       case|j  ^ 

Since  it  is  proposed  to  derive  roll  and  pitch  moments  from  the 
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control  system  by  torquing  the  controllers  about  their  vertical 
axis    (7  axis),   it  is  convenient  to  sum  torques  in  the  coordinate 
frame  of  the  case  (gu  frame).     Accordingly, 


I 


J  GU        ^*JGU 


M  I  =  ^1  -^  W,    ^^.'AhI^,,       (Eq  6.2.2) 

easel  ^TT        ^MgU  I^  GU       J  GU  ^ 


where      H  =  H     +  H      +  H 
r  g  c 


From  Appendix  D,   Equation  D.  3.  5, 


W 


I,  GU 


pcy  +  qsy 

-PS7  "^  qcY 

r  +  7 


(Eq  6.2.  3) 


The  approximation  made  in  Section  2,  6  is  that 


H 


j  GIM 


0 

J    3 

P 


0 


(Eq  6.2.4) 


Now,   in  the  gu  frame 


H 


rJ^^  =  QGU,GIMHrJ^j^ 


J    /3Cq' 
P 


(Eq  6.  2.  5) 


Define    J  as  the  combined  moment  of  inertia  of  the  case, 

cgr 
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gimbal  and  rotor. 


J  =  J      +  J      +  J 

cgr  c         g  r 


(Eq  6.  2.  6) 


Combining  the  last  two  equations  gives 


"J 


GU 


P 


p"^  cgr  ' 


(Eq  6.  2.  7) 


Differentiating  with  respect  to  time  where  it  is  assumed  that 
angular  momentum  about  the  spin  axis,   Jp^  ,    is  constant  gives 


"J 


GU 


J  l3Saa 


cgr 


(Eq  G.  2.  8) 


Then  substituting  the  appropriate  values  into  Eq  6.  2.  2  gives  for 
the  Number  1  Controller  the  following. 


^y       case  1 


•.       • 


GU 


(R   h7)J    3CaH/-pST-'-qC7)(Jp/5SQ'l=J^.^^,7) 


J  fiScva-   (pC7  +  qS7)(J   /3  Sc>+  J^.^^.  7 ) 


+  J^^^T  ^=  (^    +  pCT  ^qST)  Jhca 
*-&^  p 


(Eq  6.  2.  9) 
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The  X  and  Y  components  of  the  torques  in  the  last  equation  are 
normal  to  the  single  degree  of  freedom  of  the  case,   and  there- 
fore are  transmitted  through  the  bearing  to  the  vehicle.     The  Z 
component  is  aligned  with  the  bearing  axis,   and  the  rotation  of 
the  case  must  be  controlled  by  a  torque  generator  about  the  y 
axis.     The  single  Z  component  can  be  written  as 

••  •    • 

M  i=J  y  ■¥  J    ^  a   +  pC7   +  qST    Ca  (Eq  6.  2.  10) 

T  case  1  cgr  p^  t^^  i        ^^  i    \^u  \    m 


M  .    consists  of  a  torque  generated  part,    a  load,    and  a 

viscous  friction.     Applying  these  to  equation  6.  2.  10  gives  for 


each  of  the  controllers  the  following. 


Controller  1        Let    F,    =  J  "i  +  C 

1  cgr  dt 


^TGl    "   "^Ll    -^1  ^1    "  '^p^ 


M^^j    +   M^j    =  Fj  7i    +   Jp^^i  ia^    +  pCT^    +  qS7i)Q^i 

(Eq  6.  2.  11) 


Controller  2       {y^   =  7^    +   180°) 


%G2   +  %2   =  ^2  ^2   -^^  Jp,^2<^2  "  ^^^2  '  ^273)0^.2 


(Eq  6.  2.  12) 


Controller  3       (7     =  7,   -  90°) 


•  • 


^TG3   -^  ^LS  =   ^3  ^3   -^  \^h  <^3  +  P^^S  "  ^=1^73)0.3 


(Eq  6.  2.  13) 
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Controller  4        (7^   =   7^  +   90°) 


^TG4   "^  ^L4   =   ^4  ^4  ■*■  "^p    ^4  ^^4  '  P^'^4''  '^^^A^^'''^ 

(Eq  6.  2.  14) 

Combining  Controllers  1  and  2  in  accordance  with  the  logic  that 
7     =   -  7      and    cc-    -  a^  ,   and  making  the  further  stipulations 


that 


^L2   =   ° 


Mli    =   M^^,2 


\^-\h 


F         =       F 
1  2 


(Eq  6.  2.  15) 


which  implies  that  the  two  controllers  are  geared  together  gives 


M            -    2  J    i3p.J7iC« 
7^    =    .^^iiii E_         1  (Eq6.  2.16) 


The  assumption  that  a^    =  o-^   is  not  exact  because  the  a^  and  a^ 

gimbal  axes  cannot  be  easily  geared  together,   but  if  the  variables 
are  enclosed  in  an  error  position  control  loop  the  assumption  is 
believed  valid.     Note  that  the  gearing  of  controllers  1  and  2  com- 
pletely eliminates  the  cross  coupling  effects  of    q    and    a  ,   and 
leaves  only  the  primary  attitude  rate  variable,   p  ,   of  the  loop. 
This  is  an  important  advantage  of  a  system  using  two  controllers 
operating  back- to-back. 
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In  arriving  at  the  result  of  equation  6.  2.  16  is  must 
not  be  overlooked  that  an  important  assumption  is  that  the  angles 
a^   and  a^  are  maintained  approximately  equal.     This  assumption 

gives  the  two- degree- of- freedom  controllers  the  characteristics 
of  a  single-degree-of-freedom  controller  in  that  the  a  axis  is 
approximately  rigidly  restrained  unless  all  four  of  the  controllers 
move  in  the  same  direction  as  a  unit.     In  principle  this  may  be 
accomplished  by  wiring  the  torque  generators  controlling  the 
a  axis  in  parallel  and  requiring  the  torque  generators  to  have  a 
high  back  emf. 

Combining  controllers  3  and  4  under  the  conditions  that 
the  two  are  geared  together  and  driven  by  controller  3  gives 
for  similar  assumptions  to  equation  6.  2.  15  that 


M^^^   4=   2  J    (3qCy^Ca 
73= ^^=^ 2F^ (Eq  6.2.17) 


Note  that  in  this  equation  for  pitch  control  that  the  cross  coupling 
effects  of    roll   and  yaw  do  not  appear  in  the  equation. 

6.  3  Response  of  the  Controllers  to  Moments  Applied  by 

the  Yaw  Torque  Generators 

Following  the  procedures  of  the  previous  section 

Then  it  may  be  written  that 

ZM        I  =  hI  +  W  "^  hI 

gimj  Qjj^        ^  GIM  ^    "^I,  GIM     "J  GIM  (Eq  6.  3.  2) 


where    H  =  H     +  H 
r  J 
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Now  from  Appendix  D  the  angular  velocity  is 


W 


I.GIM 


PC7  +   qS7   +   a 


(-  pS7   +-  qC7)Cc^   +    (r  +  7)  Sc^ 


(pS7   -  qC7)Sc^   +    (r  +  y)Ca 


(Eq  6.  3.  3) 


Using  the  approximation  of  equation  6.2.4  one  would  first 


guess  that   H 


rg 


vanishes;  however,    such  is  not  the  case 


GIM 


and  we  must  use  the  following  from  Appendix  E 


Hj 


GIM 


J 

rg 

an  — 

a 

0 

0 

(Eq  6.  3.  4) 


This  step  can  be  justified  by  considering  the  complete  equations 
in  Appendix  E. 


^'"^  U  GIM 


J       a  -     ipSy  -  qC7)  S»  +  (r  +  7)Ca     J^jS 
rg  i     I        ±     .  p 


0 


[pC7   +   qS7   +   Q-]  J    ^ 


(Eq  6.  3.  5) 


The  Y  component  vanishes  only  for  the  assumptions  made  herein, 
which  means  there  are  no  large  torques  transmitted  about  the 
spin  axis  of  the  rotor.     Here  we  are  interested  only  in  the  X 
component  since  the  Z  component  will  be  transferred  through 
the  gimbal  to  the  case  where  a  component  of  it  will  appear  at  the 
7  axis.     Note  that  the  cosine  a   component  can  be  identified  in 
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the  corresponding  portion  of  equation  6.  2.  9. 

Unfortunately,  there  does  not  appear  to  be  a  simple 
method  of  mechanically  connecting  the  gimbals  of  the  four  con- 
trollers so  that  each  of  these  angles  are  exactly  the  same. 
Therefore,   each  of  the  gimbals  must  be  handled  separately. 
These  may  be  written  as  follows 


Controller  1        Let  G^    =  J^^'dt  "^  ^ 


--^TGi  "^  '^p^  ^P^^i  'q^'>'i^^i  "^  (^  +  ^1)0^1 


M. 


^1    = G^ 

(Eq  6.  3.  6) 

Controller  2 


^TG2   "^  "^p^    (pS7i  +  qCTi)Sa2    +    (r-  7^)0:^2 

^2   = G^ 

(Eq  6.  3.  7) 

Controller  3 


^TG3   '^  -^p^   ^^^'^3  '  qS'>'3^^3   +   (r  4-  73)0:^^ 

(Eq  6.  3.  8) 


-3=  G3 


Controller  4 


M 


TG4    ^  "^P^    <-pC73  -  ^373)3^4  +    (i-  -  73)Ccv 
a^ ^ 


^  (Eq  6.  3.  9) 
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In  the  above  four  equation  it  is  assumed  that   To  ~  "'^1    '  "^4  ~  ''^? 

and  that  J    S  is  held  constant. 
P 

If  it  is  assumed  momentarily  that  the  a  gimbal  axes 
can  be  mechanically  geared  together  then  one  may  sum  the  yaw 
moments  of  the  four  controllers.     Since  the  angular  momentum 
of  each  controller  is  identical  a  summation  of  yaw  moments  can 
be  accomplished  by  first  finding  the  average  gimbal  rate  by  sum- 
ming the  four  gimbal  rates.     This  gives  the  following  equation. 

» 

a= 1^ P  ^J ^ (Eq6.  3.  10) 


In  this  expression  the  moment  applied  by  the  torque  generator, 
Mrr^r^  ,   is  applied  to  each  of  the  four  gimbal  axes.     Note  that 

unlike  the  previous  equations  for  7..    and  7^  the  cross- coupling 

rates  p  and  q  do  not  vanish  although  they  are  multiplied  by  the 
product  of  two  angles  which  may  remain  small.     However,   the 
large  cosine  components  of  the  cross- coupling  rates  do  drop 
from  the  equation  as  do  the  terms  in   7..    and   7^  . 

It  is  proposed  to  gear  the  a  gimbals  together  electrically 
by  finding  the  average  gimbal  angle  and  comparing  this  with  the 
angle  of  each  gimbal  separately  in  an  error  feedback  loop.    The 
equations  for  this  operation  are  contained  in  the  next  section. 
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6.  4 Torque  Generator 

The  control  power  requirements  for  the  gyro  controllers 
are  insignificant  compared  to  the  power  required  to  drive  the 
rotors  at  operating  speed  as  will  be  shown,   and  since  the  peak 
torque  requirements  are  small,   either  a  direct  current  motor 
or  a  two  phase  alternating  current  motor  is  applicable.    The 
ideal  torque-speed  characteristics  of  both  types  are  the  same 
to  a  first  approximation. 

J^  0  +  B  0  +  M^  =  K^  V^  (Eq.  6.4.1) 

where  6      =  motor  position  angle 

J^  =  motor  moment  of  inertia 

B     =  viscous  friction  constant 

My  =  load  torque 

Krp  =  torque  constant 

V^  =  input  control  voltage 

The  motor  is  geared  to  the  controller  with  a  gear 
ratio  p  >  1;  therefore,  the  torque  applied  to  the  controller 
gimbal  is 

M^^  =  pMj^  (Eq.  6.  4.  2) 

The  angular  relation  between  the  motor  position  angle 
and  gimbal  position  angle  is 

6  =  pa     for  yaw  control 

0  =  P7     for  roll  and  pitch  control 


(Eq.  6.  4.  3) 


This  gives  for  yaw  control 


(Eq.  6.  4.  4) 
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and  for  roll  and  pitch  control 

^TG  "  ^T  P  "^M"  C^M  P^  "  +  B  p^  7}  (Eq.  6.  4.  5) 

Combining  equations  6.  2.  16  and  6.  4.  5  gives  for  roll  control 

K     pV       -  2  Jp^'pCT.  Ca 
T^  =  —^ ^    p       ^ (Eq.  6.  4.  6) 

It 


where 


^l7=0jcgr+JMp')il<2C+Bp2) 


Likewise  for  pitch  control 


.         K^pV^^2JpgqCT3Ca 

73= p (Eq.  6.  4.  7) 

37 


For  yaw  control  the  four  controllers  give  separately 
the  following  relations. 

^^'      Gl.  =  (jrg+JMP'>^<C+Bp') 
.  K^.pV^  -+ Jp^  [(pS7i-qC7i)sc.^  ^(^ +'>'l)^^l] 


^1= u 


la 

(Eq.  6.  4.  8) 


K^PV^  +Jp/3  [(pS7i  +qCyJsa^  =^(^-^l)c^2] 


""2 G. 

Za 


(Eq.  6.  4.  9) 
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•         ^TP^M^Jp^[CP^^3-^ST3)scv3+(r +73)^.3] 


0*0    = 


3  ~  G„^ 

6a 


(Eq.  6.  4.  10) 
•         ^T  P  ^M  ^  Jp  ^  [("P  ^"^3  "  ^  S'^3)s^4  +  Q  -  *4)^^4] 


Q"^  = 


40- 


(Eq.  6.  4.  11) 


In  the  above  equations  subscripts  have  been  omitted,   but  each  of 
the  equations  requires  its  own  characteristic  parameters. 

Since  all  controllers  are  slaved  to  the  position  angle  of 
the  average  controller,   the  relation  for  the  motor  voltage  of  a 
typical  controller  is  as  follows. 

VM  =  VMz-K9p(«i-a^,J  (Eq.6.4.12) 

where     ^yry  =  error  voltage  and  KO  =  feedback  constant. 

Substituting  equation  6.  4.  12  into  the  above  four  equations  gives 
the  following  when  written  in  terms  of  position  angles. 

Controller  1 

(Eq.  6.  4.  13) 

Controller  2 

^         ^T  P^MZ  ^  ^0  P  "ave  ^ ^P  ^[(P  ^^1  ^  ^  ^^i>"2<"-  ^i>"2] 

tCO^k^p 


_d 
d 


(Eq.  6.  4.  14) 
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Controller  3 


^T  0  ^MZ  +Kg  p  a^^^  +  Jpe[(p  Cy^qSy^a^^^^a^] 


(Eq.  6.  4.  15) 


Controller  4 


^TP^MZ^^eP^ave^Jp^[(-P^^3-q^^3>4<^-^3)^^4] 


dt    ^4a.^^0P 


(Eq.  6.  4.  16) 


If  we  assume  that  the  gimbal  angles  are  held  approximately- 
equal  such  that    Ca.   =  Co-^  =  Ca^  =  Co- .  ,   then  the  operation  of 
finding  the  control  moments  contributed  by  the  four  controllers 
can  be  accomplished  by  first  finding  the  average  value  of  the 
gimbal  rates  by  combining  the  last  four  equations.    This  gives 
for  the  average  gimbal  rate  the  following 

^       2K^pVj^^  +Jp/3[(pSTi-qS73^Sc^  +  2rCc^] 


2G 


(Eq.  6.  4.  17) 
Where      G^  =  (j       +  J^  p')  ^  +    (c  +  B  p^) 
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6.  5 


Attitude  Sensors 


The  primary  source  of  guidance,   navigation,   and  attitude 
information  for  a  spacecraft  is  from  external  sightings.     As  was 
shown  in  Chapter  4  the  ambient  fields  around  a  spacecraft  are 
too  feable  to  drive  a  practical  sensor,   and  this  implies  that  some 
type  of  star  tracker  must  be  provided  if  the  spacecraft  is  to  be 
self  sufficient.     A  control  system  can  control  the  attitude  of  a 
spacecraft  only  to  the  precision  provided  by  the  attitude  sensor, 
and  the  design  of  a  highly  linear  sensor  is  a  problem  within  it- 
self.    As  done  by  other  investigations  of  hardware  for  applying 
torques  to  a  spacecraft,  this  thesis  assumes  a  linear  sensor 
for  each  of  the  three  attitude  reference  axes  of  the  spacecraft. 
Accordingly  the  following  equation  is  assumed  to  represent 
the  sensor. 


K       K 
sx     ex 


0 


K       K 

sy     ey 


0 

<l>^-<l> 

V 
ex 

0 

d    -  G 
r 

= 

V 

cy 

K       K 
sz     ez 

^^-^ 

V 
cz 

— 

—            _ 

__       ^_ 

(Eq  6.  5.  1) 


K     = 
s 

K     = 
e 

<l>      = 
r 

(/,      = 
V     = 


sensitivity  of  sensor 

sensitivity  of  error  amplifier 

reference  angle 

vehicle  attitude  angle 

voltage  signal  to  tandem  compensation 


Any  practical  sensor  has  linear  characteristics  over 
only  a  small  range,    and  the  effects  of  saturation  will  be  consid- 
ered in  a  later  section  concerning  the  position  control  mode. 
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6.  6  Tandem  Compensation 

To  accommodate  the  provision  for  shaping  the  error 
signal  before  it  is  fed  to  the  control  system  torque  generators, 
a  transfer  function  which  is  presently  undefined  is  inserted  into 
each  control  loop.    The  following  equation  applies. 


X 


0 
C 


Y 


0 


0 


0 
0 

c. 


"^cx 

V 
^MX 

^cz 

= 

"^MY 

^cz 

V 
MZ 

V 


(Eq.  6.  6.  1) 
=  voltage  signal  to  tandem  compensation 
C      =  transfer  function  of  conapensation 
^M  ~  voltage  signal  to  torque  generators 


It  is  not  anticipated  that  any  tandem  compensation  will 
be  required  in  the  examples  of  this  thesis.    In  an  actual  problem 
certain  regions  of  the  s-plane  may  be  denied  to  the  control  sys- 
tem designer  because  of  unstable  modes  of  the  airframe,   etc,  , 
and  for  this  reason  space  in  the  equations  has     been  provided 
for  tandem  compensation. 
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For  convenience  we  write  equation  6.8.  1  in  shorthand  notation 


as  follows, 


2m(12-34-1234>J^  =  2h[aJ 


e    ■ 

r 

-  0 

^r 

-  i' 

+  4 


h^b] 


(Eq  6.8.2) 


where    I  A  J    and    I  BJ    are  the  indicated  matrices  and  (a..)  and 
(b..)  are  typical  elements  of  these  matrices. 

To  obtain  the  equation  for  the  complete  closed  loop 
system  we  substitute  equations  6.8.2  and  2.4.5  into  equation 
2.3.3  which  gives  the  following. 


IIxJ-2^4h\i^ 


4H^b 


12  dt  (' 


2Ha 


11 


+ 


h\ii 


h«\i^ 


(I     -  I    )  qr 

z        y'^ 

(I     -  I    )  pr 
X         z    ^ 


(I     -  I    )  qp 
y        X  ^^ 


\i'*^\2£ 


dt 


2Ha 


22 


^''\2Tt 


^'^'^isi 


-2Ha 


13 


^^^SsW 


-2Ha 


23 


/I     d_   +         2        _d_ 

z  ,^2       ^"   "^33  dt 
dt 


-2Ha 


33 


J 


~-     I^ext)^-I^IVMJ^-2H[A] 


e 

r 

^r 

(Eq  6.8.3) 
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Equation  6.8.2  represents  the  complete  loop  of  the 
attitude  control  system,   and  as  can  be  seen  there  is  non-linear 
interaxis  coupling  between  the  three  loops  as  well  as  non-linear 
terms  arising  from  the  Euler  inertia-cross  coupling.     If  it  is 
assumed  that  compensation  is  applied  to  the  loops  as  set  forth 
in  section'  3.  6,    and  also  compensation  for  the  inertia-cross 
coupling  terms,   the  following  expression  is  found  for  the  roll 
attitude  angle. 


dt 


I    ^  +  4H^b..-^  -  2Ha,,  (/)  =  M     -2Ha,,(/)  (Eq6.8.4) 

X  ,^2  11  dt  11^  X  11  ^r  ^ 


Similar  expressions  can  be  written  for    6  and  4^ . 

If  the  angular  motions  of  the    a    axes  of  the  controllers 
are  restricted  to  small  angles  then  the  following  equations  apply. 


I    F-    ^  +  (2HCtJ^-#  +  2HC7.K  (^  =  F-tM     +2'RCy,K6 
X      l7_,j.2  '  1      dt  '  1    x^  1'      X  '1    x^r 

(Eq  6.8.5) 


■^dt' 


2 

I..F..  ^-^+  (2HC7o)^  +  2HCToK   0  =  F.tM     +  2  HCToK   6 
y     iy  ^^  '3      dt  '3     y  3  '      y  '3    y   r 

(Eq  6.  8.6) 


2 

I    G    ^  +  (2  HC«)^  +  4  HCc^  K  L^  =  Ga  M     +4  HCcv  K  ^ 
z     a  .2.  dt  z^  z  z^] 


(Eq  6.8.  7) 
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where       F^^  =   (2  J^g^  +  Jm  P^>  dt    '   '^C   +  Bp^) 


^37=   <2J,g,+JMp')^^    (2C   H-Bp2) 

G       =    (J        +J„p')4'<^^^^'' 

a  gr  M  ^      dt 


the  constants  for  the  different  axes  may  not  be  equal. 

Block  diagram  representation  for  the  above  equations 
are  contained  in  Figure  6.8.1  and  Figure  6,8.2.     Figure  6.  8.  1 
contains  the  roll  and  pitch  channels  whereas  the  yaw  channel 
is  shown  in  Figure  6.8.2  in  two  parts.     Part  (a)  shows  the  four 
individual  controllers  driven  by  the  same  error  voltage.     To 
keep  these  controllers  to  the  same  angle  their  output  angles  are 
averaged  and  then  compared  with  the  actual  output  gimbal  angle. 
The  error  is  then  fed  back  to  the  gimbal  torque  motor.     Part 
(b)  of  Figure  6.8.2  is  the  result  of  ideal  controllers  and  repre- 
sents equation  6.  8.  7. 
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^$HZiHZ!lH" 


Cz 


^e-453-^4-[7^ 


^v^ 


♦  "AVE  f'EEOBACK 


HCa 


dt 


"exf 


■-4i 


H=J-/9 


%-E5He-/^ 


^^ 


"AVE  FEEDBACK 


at 


-^5. 


fO; 


«-i^;-E7p]-x4^/^ 
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o^yg   FEEDBACK 


HCa 
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ra-^ 
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fi 


—I  HCa  Uhi 
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-i® 


PART   (o)    INDIVIDUAL  CONTROLLERS  AND  AVERAGING  CIRCUIT  SHOWN 
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"ext 


f- 


HCa      ■• 


I" 


PART  (b)      SIMPLIFIED  REPRESENTATION  OF  EQUATION  6.8-7 

Figure  6.8.  2     Block  Diagram  Representation  of  the  Position 
Control  Yaw  Equation  for  Uncoupled  Motion. 
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6^ 


Rate  Control  Equations 


The  purpose  of  the  rate  control  is  to  provide  the 
astronaut  with  a  means  of  slewing  the  vehicle  at  a  particular 
attitude  rate.     Since  the  position  of  the  gimbal  angles  are 
directly  related  to  the  amount  of  control  system  angular  mom- 
entum transferred  to  the  vehicle,   the  gimbal  angles  represent 
attitude  rates.     Therefore  the  rate  control  is  essentially  an 
error  feedback  control  loop  to  control  the  gimbal  angle,    and 
the  block  diagram  of  the  roll  channel  is  shown  in  Figure  6.  9.  1. 
During  this  operation  the  error  signals  from  the  vehicle  attitude 
sensors  are  open  circuited. 

During  the  rate  control  maneuver  the  voltages  to  the 
torque  motors  are  given  by  the  following  equation. 


V 


1    !" 


mx 


V 


my 


V 


mz 


K 


Tl 


0 


0 
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Ir 
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7c 
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L. 


(Eq  6.  9.  1) 


The  voltages  given  by  equation  6.9.1  may  be  substituted  into 
equations  6.4.6,    6.4.7,   and  6.4.  17  to  find  an  expression  for 
the  gimbal  angles  as  a  function  of  the  rate  variables  and  the 
rate  inputs.     For  shortness  we  write 

1  7I     tx^x 


2  73     ty^^y 


K,    -   K    K^    p 
o  en     ten  "^a 


which  gives  the  following  matrix  equation. 
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1 

_. 

l7dt             1 

0 

0 

^1 

0 

^37^    "   ^ 

0 

^3 

0 

0 

^dl  ^  s 

Q' 

1  '  Ir 


^2^3r 


^3-r 


+  2H 


-CTjCcc 


S7-  Sa 


+  CT3CQ'  0 


S70  Sa  Co- 

— ^ +  — 


l_  -» 


(Eq6.  9.2) 


The  above  matrix  equation  must  be  solved  for  the  gimbal  angle 
rates  so  that  these  rates  can  be  applied  to  the  control  system 
represented  by  equation  G.  2.  11.     Because  of  the  non-linearities 
caused  by  the  sine  and  cosine  terms  of  equation  6.  9.2,   the 
Laplace  techniques  are  not  simple  to  apply  unless  the  angles  can 
be  considered  constant.     Therefore,   in  the  equations  that  follow 
the  notation  employs  the  letter  s  so  that  a  closed  form  solution 
can  be  written,   but  this  symbol  s  refers  to  the  Laplacian  oper- 
ator only  in  the  special  case  where  sine  and  cosine  values  of 
the  gimbal  angles  can  be  considered  constant.     Otherwise,   for 
cases  where  the  gimbal  angles  make  large  changes,   the  symbo] 
s  must  be  interpreted  as  a  differentiation  with  respect  to  time. 

Solving  the  matrix  equation  6.  9.  2  for  the  gimbal  angle 
rates  and  substituting  in  equation  G.  2.  11  gives  the  following 
expression  for  the  moments  generated  by  the  spacecraft  control 
system. 
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As  in  the  previous  solution,    rewrite  equation  6.  9.  3  as  follows, 


SF\     +lc7)'^lr 

I7  1 


ZCaCjo  SKp 


Control!  A 


37      "2 
2  0^  SK. 


V^ 


SG     +  K. 


o: 


4  H^  [e] 


(Eq  6.9.4) 


Substituting  equation  6.9.3  and  2.4.5  into  equation  2.  3.  3  gives 
the  following  equation  for  the  response  to  rate  control. 
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(Eq  6.  9.5) 
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The  uncoupled  set  of  equations  for  roll,   pitch  and  yaw 
is  given  by  the  following. 


S(I    F,     S    +  I   K      +    (2HC7i)^] 
X      I7  X     1  1      / 


(Eq  6.  9.  6) 


(Eq  6.  9.  7) 


M    (SG      +  K^)    -  ai^HCaK^S) 
z a o J 

S(I  G    S   +  I   K„    +   (2HCTo)^) 


(Eq  6.  9.8) 


The  following  equations  are  valid  for  the  rate  input 
commands  which  drive  the  gimbal  angles  to  large  values  in 
those  cases  where  there  is  negligible  inter-axial  coupling. 

For  Roll 


■N 


F,^,  T,  +  K,7,    =   K,7,^   -    2HC7iP 


I7   '1  I'l 


I'lr 


M      =    -2HC7,7,    +  I    p 
X  '11  X  ^ 


(Eq  6.9.  9) 


For  Pitch 


>, 


^37^3   ^   ^-^3   =   ^2^3r   "^  ^  HC73  q 


M 


y 


2HC7373    +1    q 


(Eq  6,  9.  10) 


For  Yaw 


Go-    +   Kocc 


S-r  ^  «^r 


M     =   4HCaa   +  I    r 
z  z 


(Eq  6.  9.  11) 
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For  the  completely  general  case  it  is  necessary  to 
solve  equations  6.  9.  2,   G.  2.  11,   2.4.  5,   and  2.  3.  3  by  machine 
computation. 
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Figure  6.9.1    Block  Diagram  Representation  for  Roll  Channel 
of  the  Uncoupled  Rate  Control  Equation, 
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6.10  Summary 

This  chapter  has  served  to  derive  the  closed  loop 
equations  of  motion  of  the  spacecraft.     Examination  of  the 
closed  loop  position  control  equation  6.  8.  3  and  the  closed 
loop  rate  control  equation  6.  9.  5  confirms  the  fact  that  the 
equations  are  coupled  from  control  interaction  as  well  as  the 
Euler  rate-product  terms.     Therefore,   a  completely  general 
solution  requires  assistance  from  a  computer.     One  of  the 
reasons  for  choosing  the  four  controller  system  is  that  many 
of  the  large  cross- coupling  terms  in  the    primary  control 
matrix  vanish,    and  since  the  four  controller  system  is  a  zero 
momentum  type  system  the  cross- coupling  terms  of  the  gyro- 
scopic coupling  matrix  are  small.     Cross- coupling  is  a  problem 
in  any  type  of  flying  machine  if  we  expect  to    actuate  the  con- 
trol of  two  axes  to  their  full  travel  simultaneously.     Generally 
speaking  vehicles  are  not  operated  in  this  manner,    and  if 
crossed  controls  are  required  one  of  the  axes  will  predominate. 

In  the  following  chapter  the  equations  of  motion  are 
solved  to  determine  the  response  of  the  vehicle.     This  is  ac- 
complished by  assuming  no  interaxial  coupling,    and  solving 
the  equations  analytically.     Making  this  assumption,   many  of 
the  equations  can  be  solved  in  closed  form  whereas  others 
require  graphical  methods.     It  is  believed  that  this  approach 
gives  a  better  understanding  and  presentation  of  the  system 
than  a  pure  computer  solution  which  gives  graphical  response 
of  the  vehicle  to  input  disturbances. 
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CHAPTER  7 
DYNAMICS  OF  THE  SPACECRAFT 

7.  1  Introduction 

This  chapter  evaluates  the  response  of  the  spacecraft 
using  the  equations  developed  in  the  preceding  pages.    This  is  ac- 
complished by  assuming  that  the  interaxial  coupling  among  the 
three  control  axes  can  be  neglected.    This  assumption  is  valid 
provided  the  control  system  is  not  appreciably  saturated,   and 
that  the  astronauts  do  not  attempt  to  simultaneously  actuate  the 
controls  to  yield  large  rates  about  two  axes.    In  any  event,   a 
necessary  condition  is  that  the  attitude  control  system  must  per- 
form satisfactorily  about  each  of  the  individual  control  axes  so 
that  initial  analysis  of  the  equations  of  motion  on  an  uncoupled 
basis  is  considered  reasonable.    Also,    several  parameters  of  the 
system  must  be  numerically  chosen,   and  their  selection  is  greatly 
simplified  using  uncoupled  equations.    These  parameters  include 
the  damping  coefficient  of  the  controllers,   the  gain  of  the  rate 
control  loop,   and  the  gain  of  the  position  control  loop.    It  is  visu- 
alized that  the  spacecraft  will  be  operated  in  one  of  three  modes. 

1.  Zero  Input  Mode 

2.  Rate  Control  Mode 

3.  Position  Control  Mode 

In  addition  to  these  three  modes,   a  mode  called  the  Adaptive 
Mode  will  be  discussed. 

The  zero  input  mode  is  that  which  results  when  no  signal 
is  provided  to  the  torque  motors  of  the  controllers.    The  control 
system  then  operates  as  a  rate  stabilizer,   and  the  gyro  controllers 
provide  inherent  sensing.    The  stability  of  the  zero  input  mode 
will  serve  as  a  criterion  for  the  selection  of  the  damping  coefficient 
of  the  controllers.    The  damping  coefficient  is  then  held  constant 
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for  the  other  two  modes.    The  purpose  of  the  rate  control  is  to 
command  the  spacecraft  to  perform  a  particular  attitude  rate. 
The  stability  of  the  rate  control  mode  will  be  used  to  ascertain 
the  gain  of  the  rate  control  loop,   and  the  stability  of  the  position 
control  mode  will  determine  the  gain  of  this  loop.  The  adaptive 
mode  is  the  vehicle  response  to  a  controller  failure. 

The  damping  chosen  for  the  controllers  differs  appreciably 
from  that  used  in  systems  which  overdamp  the  gyro  to  the  point 
that  it  can  be  approximated  as  a  first  order  system.    In  the  present 
design  the  controllers  act  as  a  true  second  order  system  with 
two  poles  having  both  real  and  imaginary  parts.    This  of  course 
complicates  the  solution  of  the  equations,    some  of  which  cannot 
be  solved  in  closed  form.    Therefore,    in  many  of  the  equations 
which  follow  it  is  necessary  to  approximate  the  second  order 
system  by  a  first  order  system  to  arrive  at  useful  results. 

As  was  done  in  Chapter  5,   the  analysis  will  be  accom- 
plished using  the  equations  for  roll.    The  equations  for  pitch  con- 
trol are  almost  precisely  the  same  as  those  for  roll,   and  generally 
the  same  parameters  are  required  except  that  the  numbers  will 
differ  because  of  different  moments  of  inertia  of  the  spacecraft 
about  the  roll  and  pitch  axes.    The  yaw  equations  differ  slightly 
from  those  of  roll,   and  a  sumnnary  of  the  differences  will  be 
made. 
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7.  2 Zero  Input  Mode 

The  uncoupled  equations  for  the  zero  input  motion  are 

given  by  equations  6.  8.  5,    6.  8.  6,    and  6.  8.  7  where  the  position 

input  terms  have  been  dropped,   and  these  equations  are  written 

as  follows. 

F         M 

p= ^-^ ^ 5-  (Eq.    7.2.1) 

I^Fj^s  +4H^  (CT^)'' 


^37    ^v 
-^ ^-^ ^ 5-  (Eq.    7.2.2) 

IyF3^S+4H^(C73)^ 


G         M 

r  = 1 ^  (Eq.    7.  2.  3) 

I    G    s  +  4H      {Ca) 
z     a 


The  stability  of  the  spacecraft  for  the  complete  range  of 
gimbal  angles  can  be  investigated  by  assuming  small  perturbances 
wherein  the  gimbal  angle  can  be  taken  as  a  constant.    Consider  the 
equation  for  roll  motion  and  let  F^      which  is  actually  equal  to 

^1     "^    C^*^      r  *  ^M^O^  +  (2C  +Bp^)    be  written,    for  short- 
ness,   simply  as      F^     =  J  s  +  k.     This  gives 

M  (s  +  k/J) 

p= 7y 5 (Eq.    7.2.4) 

2       ks        4H^(C7i)^ 

I         s^+^+ r— F 

X  J  I    J 

X 
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The  motion  follows  that  of  a  second  order  system  with  a  damping 
ratio 


I  = 


4HCTi    V       J 


I 


(Eq.    7.  2.  5) 


and  a  natural  frequency. 


2HCt 


1 


o_  = 


X 


(Eq.    7.  2.  5) 


The  effect  of  the  disturbance  is  least  when  the  gimbal  angle  is 
small  for  which  the  following  are  defined 


^o       4H    V     J 


(Eq.    7.  2.  6) 


GO        = 

o 


2H 


X 


(Eq.    7.  2.  7) 


The  solution  for  the  response  of  the  system  to  an  impulsive  mo- 
ment,   M   ,   for  tl 
o 

lowing  equation. 


ment,    M   ,   for  the  lightly  damped  controller  is  given  by  the  fol 


M 


o_ 

X 


ct. 


,L    2    „2  ^ 


e  .    sin< 


^o  ^CTi-C^t^  +  tan  "  ^ 


^. 


Ct?-^^ 


O 


(Eq.    7.  2.  8) 


Examination  of  this  equation  shows  that  the  roll  rate  is  always 
stable  for  any  gimbal  angle  less  than  90  degrees  because  the  ex- 
ponential damping  factor  is  not  a  function  of  the  gimbal  angle. 
Critical  damping  occurs  when    CTh   =  C       and  the  solution  for  this 
motion  is  given  by  the  following  equation. 


176 


Mr  -       -?;    oj    t 

O  \     y  4.  ,  1  O  O 

p  =  -f —   K    CO    t  +   1       e 
^        I        L^o    o  J 


(Eq.    7.  2.  8) 


The  highly  damped  motion  can  be  approximated  by  assuming  that 
the  gimbal  moment  of  inertia  is  negligibly  small  giving  the  fol- 
lowing equation  for  roll  rate. 

M                        1 
x_     

.       .  (Eq.    7.  2.  9) 

The  solution  to  this  equation  for  an  impulsive  moment  disturbance 
is  given  by 

M  -^h'ct^^ 

O 


p  =^ —    e         Ik 

^x  ^ 

^  (Eq.    7.  2.  10) 

In  terms  of  previously  defined  parameters  this  motion  can  be 
written  as  follows. 


""       t 


M 


21, 


p  =  -j^    e  (Eq.    7.  2.  11) 


X 


The  analysis  of  the  response  due  to  a  step  moment  dis- 
turbance is  much  more  difficult  to  analyze  because  the  gimbal 
angle  cannot  be  assumed  constant.   As  a  beginning,    however,    let 
us  assume  that  the  moment  is  small  enough  that  first  order  dyna- 
mics can  be  used  in  which  case  equation  7.  2.  9  is  considered. 
Further,    let  the  gimbal  angle  remain  small  so  that  the  cosine  of 
the  gimbal  angle  can  be  approximated  as  a  value  of  one.    A  step 
of  M.  at  time  zero  gives  the  following. 
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M,  k 
p  =  — H-\l  -  e  I  (Eq.    7.  2.  12) 

The  final  value  is  seen  to  be  independent  of  the  moment  of  inertia 
of  the  vehicle.    The  roll  rate  of  a  vehicle  without  an  attitude  control 
system  in  response  to  a  step  disturbance  moment  is  given  by 

Ml 
p  =-p^  t  (Eq.    7.  2.  13) 

X 

The  roll  rate  of  a  spacecraft  without  a  control  system  connpared 
to  that  with  gyro  controllers  both  of  which  experience  a  step  mo- 
ment disturbance  is  given  by  dividing  equation  7.  2.  13  by  7.  2.  12 
as  follows. 


^without  control  system    _     ^o 
^with  gyro  controllers  o 


(Eq.    7.  2.  14) 


''^o 


The  factor   — —    is  the  reciprocal  of  the  time  constant  for  the 

o 
gyro  controller  and  the  factor  is  large  for  fast  controllers,    so 

that  the  ratio  is  plotted  initially  as  a  steep  curve  in  Figure  7.  2.1. 

The  effect  of  saturation  of  the  gyro  controllers  is  sketched  in  to 

show  that  the  ratio  approaches  unity  as  time  increases. 
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Figure  7.  2.  1         A  plot  of  the  ratio  of  roll  rate  of  an  uncontrolled 
spacecraft  compared  to  a  gyro  controlled  spacecraft  both  of 
which  experience  a  step  moment  disturbance.    Dotted  curve  is 
undocumented. 
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The  solution  for  a  step  moment  disturbance  wherein  the 
controller  gimbal  angles  move  through  their  complete  range  is 
solved  by  a  phase  plane  analysis,   however,   first  consider  the 
special  case  of  recovery  from  a  large  roll  rate  in  the  absence 
of  external  moments.    The  solution  for  the  gimbal  angle  has  been 
determined  as  follows. 


00 


2C 


7.  =  arc  tan  < 


tan  7 


(Eq.  7.  2.15) 


Using  this  relation  for  the  gimbal  angle  gives  the  following  equa- 
tion for  the  roll  rate. 


p  =  p  sin 

'^      ^max 


arc  tan^  e 


2;, 


tan  7 


o 


J(Eq.  7.  2.  16) 


These  equations  are  plotted  in  Figure  7.  2.  3  for  an  initial  gimbal 
angle  of  89  degrees. 


r,    600- 


Figure  7.  2.  3  Plot  of  gimbal  angle,    7.,   and  non-dimensional 

roll  rate,   p/p  ,    for  spacecraft  controlled  by  two  gyro  control- 

lers  with  initial  condition  of  7,(0)  =  89  degrees.    Total  angular 
nnomentum  of  system  is  zero  and  there  are  no  external  moments. 
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This  special  case  of  recovery  from  a  large  roll  rate  in 
the  absence  of  external  moments  also  has  a  simple  phase  plane 
solution  as  follows.   As  shown  in  Figure  7.  2.  4,   plot  the  following 
equations  for  selected  values  of  gimbal  angle. 

.        (2HC7i)^ 

p  +  ^ j^-*—    p  =  0  (Eq.  7.2.17) 


X 


ffi'     'p  =  p  STi  (Eq.  7.  2.  18) 

^         ^       ^max       1  ^ 

The  trajectory  is  determined  by  intersections  of  lines  of  the  same 

gimbal  angle.    The  times  of  the  trajectory  can  be  found  from  the 

ratio      Ap/p  ,   and  from  this  the  curve  roll  rate  versus 

^  ^average 

time  can  be  plotted.    Thereafter,   roll  rate  can  be  integrated  to 
find  the  roll  displacement  as  follows. 

0  =     \  pdt  (Eq.  7.  2.  19) 
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Figure  7.2,5      The  relationship  between  gimbal  angle  and  roll 
displacement  for  the  zero  input  mode  with  zero  initial  conditions 
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A  closed  form  solution  of  equation  7.  2.  1  is  difficult  even 
when  only  the  first  order  dynamics  are  considered.    There  is, 
however,   a  direct  relationship  between  the  vehicle  attitude  angu- 
lar displacement  and  the  gimbal  angle.    For  example,   the  follow- 
ing equation  is  always  true  for  zero  initial  conditions  and  for  zero 
input  to  the  controller. 

JT*^  +  kT=  -2HC7P  (Eq.  7.  2.  20) 

Considering  only  first  order  dynamics  the  equation  can  be  inte- 
grated to  give  the  following  result. 


2_H± 


7=2  arc  tan^  e 


This  equation  is  plotted  in  Figure  7.  2.  5. 


y  (Eq.  7.  2.  21) 


An  approximate  phase  plane  analysis  can  be  made  of  the 
position  input  mode  operating  to  minimize  a  step  torque  distur- 
bance by  assuming  that   7  is  constant.    The  first  order  differential 
equation  approximation  of  equation  7.  2.  1  is  given  by 


^        4H^(C7i)^  ^v 

p+ 1     ]^  P   =  ""i  (Eq.  7.  2.  22) 

x  X 


This  equation  is  plotted  in  Figure  7.  2.  6  for  various  contours  of 
7.  which  also  represent  contours  of  constant  time  since 

A7i 
At  =     *  (Eq.  7.  2.  23) 

Since  increased  gimbal  angles  demand  increased  values  of  the 
roll  rate,   p,    it  is  clear  that  p  must  also  increase.    At  the  end  of 
the  transient  when  p  =  M^    k/4H     ,   p   is  very  small,   but  its 
average  value  over  the  interval  from  7.   =0°    to    7.   =  -  10°    can 
be  found  by  reading  from  Figure  7.  2.  6  the  value  of 


4h2 
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and  computing  p  from  the  relation 

Pave  =  if-  (Eq.  7.2.24) 

Figure  7.  2.  7  is  an  analog  simulation  of  the  zero  input  mode  and 
shows  that   7  holds  constant  at  -M/2H  until  the  gimbal  angle  gets 
beyond  about  40  degrees.    A  phase  trajectory  determined  by  the 
method  outlined  above  is  shown  in  Figure  7.  2.  6. 
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Figure  7.2.6    Estimated  Phase  Trajectory  for  Zero  Input 

Mode  Response  to  a  Step  Moment  Disturbance. 
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All  of  the  equations  for  zero  input  control  have  been  presented 
in  parametric  form  with  no  actual  numerical  values  substituted  in 
the  equations.    Therefore,   the  equations  are  valid  for  any  general 
twin  gyro  controller  system  where  interaxial  coupling  is  negligible. 
The  zero  input  mode  is  purely  a  rate  stabilization  system.    It  has 
been  shown  that  an  impulsive  torque  disturbance  acting  on  the 
vehicle  will  cause  the  spacecraft  to  roll  through  a  displacement 
angle,   but  that  the  roll  rate  will  damp  to  zero.    A  step  moment 
disturbance  will  cause  a  steady  state  roll  and  consequently  an  in- 
creasing roll  displacement.    The  zero  input  system  is  primarily 
a  rate  stabilization  system  and  the  parameter  which  can  conven- 
iently be  chosen  to  give  the  desired  response  is  the  damping 
coefficient,   k,    of  the  gyro  controller.    It  is  not  necessary  at  this 
point  to  choose  a  particular  value  of  k  other  than  to  say  that  this 
quantity  is  chosen  on  the  basis  of  equation  7.2.6  which  solved  for 
the  damping  coefficient  is  written  as  follows. 

k  =  4HC      ^f4-  (Eq.  7.  2.  27) 

o  V  I 

X 

From  a  practical  standpoint  all  of  the  quantities  in  equation 
7.  2.  22  except  t,     are  dictated  in  the  design  by  other  considerations. 
I     concerns  the  overall  vehicle  design  and  is  determined  by  the 
mission  and  the  detail  design  of  the  spacecraft.    The  angular  mo- 
mentum of  the  controller,    H,    is  determined  from  controllability 
requirements,   and  the  combined  inertia  of  the  gimbal,    case,   and 
rotor,   J,    is  to  be  determined  by  optimum  design  of  the  controller 
to  minimize  the  ratio  of  the  total  mass  of  the  controller  compared 
with  the  angular  momentum  of  the  controller.    This  leaves  the 
parameter  I,      which  is  to  be  chosen  by  the  desired  location  of  the 
closed  loop  poles  of  the  zero  input  mode.    Although  the  damping 
coefficient,   k,    is  chosen  primarily  on  the  basis  of  the  desired 
response  to  disturbances  of  the  zero  input  mode,   the  choice  does 
affect  the  other  modes.    For  example,   the  closed  loop  poles  of 
the  zero  input  mode  will  be  the  open  loop  poles  of  the  position 
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control  mode.    Therefore,   let  us  delay  the  numerical  choice  of 
the  damping  ratio,   ^^ 
have  been  examined. 


the  damping  ratio,   t,    ,   until  the  equations  of  the  other  modes 
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Figure  7.2,7    Spacecraft  Response  in  Roll  to  a  Step 
Moment  Disturbance  oi  10  Ibrft.      I     =   10^     Ib-ft-sec^, 
J  =   10  Ib-ft-sec^    H    =10"*    lb.-ft-sec,  ^^  =  0.866.     Data  from 
Analog  SXimulation. 
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7.  3       Rate  Control  Mode 

The  equations  for  rate  input  control  given  by  Equation  6.  9.  9 
can  be  solved  graphically  by  a  plot  similar  to  7.  2.  4  where  the 
trajectory  is  defined  by  the  intersections  of  lines  of  equal  ginnbal 
angle.    For  this  solution  equation  6.  9.  9  is  written  in  a  slightly 
different  fornn  which,    neglecting  the  external  monnent  disturbance, 
is  as  follows. 

■   ^   (2HC7)^        ^^  >)  ^1    2HC7 

P         I     k  P      Vlr   ^1/  I,    k 


(Eq.    7.  3.  1) 


This  equation  is  plotted  for   p  versus  p   in  Figure  7.  3.  1  and  gives 
contours  of  constant  gimbal  angle  which  slope  from  left  to  right. 
The  second  equation  needed  is  that  of  equation  7.  2.  18  which  is 
the  requirement  that  angular  momentum  is  conserved.    This 
equation  is  plotted  in  Figure  7.  3.  1  as  vertical  contours  of  con- 
stant gimbal  angle.    The  intersections  provide  the  solution.    The 
gain  of  the  gimbal  position  control  loop,    K. ,   has  been  allowed 
to  have  a  relation  to  other  parameters  of  the  vehicle  as  follows. 

K^  =   m4^  ^gq     ^    3    2) 

X 

The  factor  m  miay  be  called  a  gain  ratio,   and  Figure  7.  3.  1  indi- 
cates that  for  high  gain  ratios  the  roll  rate  is  an  exponential 
represented  by  a  constant  slope  on  the  plot  of  p  versus  p.    The 
effect  of  the  gain  ratio  factor,    m,    is  shown  more  clearly  in 
Figure  7.3.2  which  is  a  plot  of  the  steady  state  gimbal  angle 
versus  the  value  of  the  step  input.    This  plot  shows  that  for  very 
high  gain  the  gimbal  angle  is  linear  with  the  input  angle,   but  for 
a  gain  ratio  of  one  the  variation  is  highly  non-linear.    If  the  gain 
ratio  is  less  than  one,   the  position  control  system  does  not  have 
the  capability  of  holding  the  conmplete  range  of  gimbal  angles  in 
a  stable  manner.    Figures  7.  3.  3  and  7.  3.  4  further  illustrate  the 
effect  of  low  gain  ratios.    Figure  7.  3.  3  shows  a  plot  of  steady 
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Figure    7,3.1      Phase  trajectory  for  spacecraft  response 
to  a  full  step  displacement  rate  control  input.     The  external 
moment  disturbance  is  zero.     The  response  is  defined  by 
the  intersections  of  the  lines  of  constant  gimbal  angles 
plotted  for  the  two  equations  of  motion  shown  on  the 
figure. 
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Plgure  7,3.2      Plot  of  steady  state  gimbal  angle,   7p 
resulting  from  step  input,   7,    ,  for  rate  input  mode. 
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Figure  7.3,3      Steady  state  values  of  roll  gimbal  angle  for 
gain  ratios  less  than  one  in  the  gimbal  angle  position  control 
loop  used  in  rate  input  control. 
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Figure  7,3.4  Plot  of  roll  capability  for  gain  ratios  less 
than  one  in  the  gimbal  angle  position  control  loop  used  in 
rate  input  control  mode . 
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state  gimbal  angle  for  a  step  input  of  ninety  degrees  versus  the  full 
range  of  gain  ratios  less  than  one.    Careful  inspection  of  the  cur- 
vature shows  an  inflection  point  at  45  degrees  which  was  shown 
to  be  a  peak  control  power  point  in  section  5.  31.    Figure  7.  3.  4 
shows  the  roll  capability  for  a  step  input  of  ninety  degrees  for 
the  range  of  gain  ratios  less  than  one. 

From  a  stability  standpoint  the  rate  input  mode  has  a  natural 
frequency  from  equation  6.  9.  6  equal  to 


CO 

n 


(Eq.    7.  3.  3) 


This  natural  frequency  is  greater  than  that  for  zero  input  control 
given  by  equation  7.  2.  5. 

The  product  of  the  damping  ratio  and  the  natural  frequency 
are  the  same  for  zero  input  and  rate  input.    Thus 

2^^  00^   =  -J-  (Eq.    7.3.4) 

If  the  damping  coefficient,   k,    is  chosen  on  the  basis  of  zero  input 
control  then  according  to  equation  7.  2.  6, 

X 

Substituting  this  k  into  equation  7.  3.  4  gives 

4H^ 
2C    oo„  =     I ^^  (Eq.    7.  3.  6) 


r 


"^/77 


Then  further  substitution  of  K     in  equation  7.  3.  2  into  equation 
7.  3.  3  gives 


c^^=     J    ^"^^  Vm  +(Ct)^"  (Eq.    7.3.7) 


xJ 
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Therefore, 


^  ^     (Eq.    7.3.8) 

r         I TV 


I ,^   .2 

V  m  + 


(Ct) 

which  is  to  say  that  the  damping  ratio  for  rate  control  is  propor- 
tional to  the  damping  ratio  for  zero  input  control  by  the  factor 

1  /V  m  +  (C7)       .    The  greatest  difference  between  the  two  damping 
ratios  occurs  at  zero  gimbal  angle.    Hence, 

^^ 
C^  ="1 V  (Eq.    7.3.9) 

^       Vm  +  1 

Actually,    large  values  of  the  gain  factor  m  gives  more  precise 
positioning  of  the  gimbal  angle  as  shown  by  Figure  7.  3.  2,   but  the 
greater  this  gain  the  more  lightly  damped  (more  oscillatory)  be- 
comes the  rate  control  mode.    It  appears  that  an  acceptable  solu- 
tion lies  in  the  choice  of  the  gain  ratio,    m,   as  unity.    For  a  gain 
ratio  of  unity  the  gimbal  angle  can  be  positioned  at  the  full  ninety 
degree  position,  yet  for  the  range  of  45  degrees  or  less  the  curve 
for  m  =  1  of  Figure  7.3.2  is  approximately  linear.    A  gimbal 
angle  of  45  degrees  gives  more  than  70  of  the  maximum  roll  rate; 
therefore,    it  is  not  likely  that  the  vehicle  will  be  operated  at 
angles  of  more  than  45  degrees  except  for  the  full  roll  rate  posi- 
tion of  90  degrees.    For  unity  gain  ratio 

^     =     A  (Eq.    7.  3.  10) 


^       V2 

By  cross-plotting  the  steady  state  response  to  an  impulsive 
disturbance  torque  on  the  spacecraft  shown  in  Figure  7.3.5  the 
roll  rates  can  be  compared  with  those  for  zero  input  mode  and 
those  for  no  control  system  aboard  the  vehicle.    For  unity  gain 
the  rates  are  approximately  one-half  that  of  a  vehicle  without  a 
control  system. 
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Figure  7.3.5     Plot  of  steady  state  roll  rate  versus 
magnitude  of  torque  disturbance  impulse  for  rate  control 
mode.     Infinite  gain  ratio  corresponds  to  controllers  that 
are  held  rigidly  fixed  in  the  vehicle,   and  is  the  same  as 
vehicle  response  with  out  a  control  system.     Zero  gain 
ratio  corresponds  to  zero  input  mode  of  section  7.2. 
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7.  4  Position  Control  Mode 

The  position  control  mode  is  formed  by  closing  an  attitude 
angle  feedback  loop  around  the  zero  input  mode.    A  convenient 
method  for  studying  the  position  control  mode  is  a  locus  of  roots 
which  is  a  useful  technique  found  in  most  control  theory  texts. 
Using  root  locus  theory  it  is  known  that  the  closed  loop  poles  of 
the  zero  input  mode  form  the  open  loop  poles  of  the  position  con- 
trol mode.    Therefore,   the  open  loop  equation  for  the  position 
control  mode  can  be  written  for  roll  motion  as 


(K    2HCTi^  4>     +  f, 
V  X ly     r V 

loop 


M 
'p  K^x^lT^^^"    (^^1)) 


where  the  gimbal  angle  must  be  considered  approximately  con- 
stant.   The  closed  loop  equation  is  given  by 


(k     2H  CTi"^  (/)     +  F.      M 
,  _  V    X 'ly    r  I7      X 


(i^Fj^s2+4h2(cTi)'s+K^2HCTi) 

(Eq.    7.  4,  2) 

The  poles  of  the  denominator  of  equation  7.  4.  1  are  to  be  deter- 
mined so  that  the  position  control  nnode  has  good  response  to  input 
commands  and  to  minimize  the  disturbances  of  external  torques 
on  the  spacecraft.   An  interesting  theorem  which  applies  to  a  con- 
trol system  with  no  zeros  and  three  open  loop  poles,    one  of  which 
is  at  the  origin,    is  as  follows:     the  natural  frequency  of  a  unity 
feedback  control  system  having  no  zeros  and  three  open  loop 
poles,    one  of  which  is  at  the  origin,    can  never  be  greater  than  the 
natural  frequency  of  the  system  operating  open  loop.    This  theorem 
is  easily  proven  by  substituting  jto    in  the  characteristic  equation 
of  the  closed  loop  system  and  equating  the  imaginary  part  to  zero. 
Thetheoremis  illustrated  in  Figure  7.  4.  1  by  noting  that  all  stable 
poles  of  the  closed  loop  system  lie  within  the  circle  of  radius  00 
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Figure   7.4.1 
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Figure    7,4.1      Illustration  of  four  possible  geometric  patterns 
for  the  locus  of  roots  for  the  position  control  mode.     The  three 
poles  shown  are  the  poles  of  the  closed  loop  zero  input  mode 
which  constitute  the  open  loop  poles  of  the  position  control 
mode.     The  closed  loop  poles  of  the  position  control  mode  will 
lie  on  the  root  loci  depending  on  the  value  of  the  gain  of  the 
position  control  loop. 
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which  is  the  natural  frequency  of  the  open  loop  system.    The  pri- 
mary purpose  in  Figure  7.  4.  1  is  to  illustrate  the  four  possible 
geometric  patterns  that  the  locus  of  roots  of  the  position  control 
system  can  have.    If  the  damping  ratio  of  the  zero  input  mode,  C    , 
is  greater  than  one,   then  Figure  7.  4.  1  (a)  applies.    If  the  damping 
ratio  is  less  than  one  but  greater  than  0.  866  then  Figure  7.  4.  1  (b) 
applies.    If  less  than  0.  866  then  Figure  7.  4.  1  (c)  represents  the 
locus  of  roots.    The  case  where  t,     is  exactly  equal  to  0.  866  is 
shown  in  Figure  7.  4.  1  (d).    The  overdamped  case  (a)  has  poor 
response  at  low  gains,   and  it  becomes  too  oscillatory  at  the  high- 
er gains.    Values  of  I,     less  than  0.  7  are  not  satisfactory  because 
the  rate  control  mode  and  the  position  control  mode  both  become 
too  oscillatory.    Values  of  t,     between  0.  7  and  1.  0  give  an  accept- 
able solution  so  that  the  vehicle  has  good  respons'e  to  input  com- 
mands as  well  as  good  response  to  damp  torque  disturbances. 
As  an  example  of  this  section  let  us  solve  the  equations  of  motion 
using  a  damping  ratio  of  ^     =0.  866  such  that  the  locus  of  roots 
of  Figure  7.  4.  1  (d)  applies.    This  is  not  suggested  as  an  optimum 
choice  but  is  considered  to  be  an  example  of  a  satisfactory  loca- 
tion of  the  closed  loop  poles.    Choosing  a  value  of  K     such  that  all 
three  closed  loop  poles  coalesce  gives  the  following  equations. 

=   0.  77    -± 

X 


K     =   0.  77    -^zz  (Eq.  7.  4.  3) 


I    J 

X 


0  = 


0.  192  Go^"^  CTi  (/)^  +  (s  +  3  (O.  577  oj  ")  )  ^x^\ 


S^  +  3(0.  577oo^^S^+  3(6.  577oo^^^T^)^S+0.  192aj^CT^ 

(Eq.  7.4.4) 

where    u^   is  natural  frequency  of  zero  input  mode  given 
by  equation  7.  2.  7. 
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To  further  explain  equation  7.  4.  4,    if  the  cubic  denominator 
is  written  in  the  following  form 


(s'+  2^pU.^S4-.^2^   (s  +  0=   0 


(Eq.    7.  4.  5) 


then  the  solution  chosen  as  an  example  has  the  following  values 
for  the  parameters. 


'p-' 


00      =  0.  577  oj 
n  o 


?     =   0.  577 


o 


(Eq.  7.  4.  6) 
(Eq.  7.  4.  7) 
(Eq.    7.  4.  8) 


The  primary  purpose  of  the  position  control  mode  is  to  enable  the 
spacecraft  to  track  a  reference  line.    This  reference  line  may  be 
the  line  of  sight  to  a  star,   and  usually  in  this  mode  the  vehicle 
attitude  rates  are  either  small  or  are  constant.    It  is  considered 
that  equation  7.  4.  4  can  be  evaluated  for  constant  gimbal  angles 
if  we  assume  that  the  error  signal  from  the  sensor  is  limited  to 
some  maximum  value.    In  actual  practice  this  limiting  will  be 
accomplished  by  the  saturation  of  the  sensor.    If  the  maximum 
error  signal  from  the  sensor  is  </>    then  the  steady  state  roll  rate 
resulting  from  the  constant  sensor  error  can  be  derived  as  the 
following  equation. 

1/2 
1-/^. 
V^m. 


^     =5.2Vf 

X 


max 


lax, 


(Eq.    7.  4.  9) 

This  curve  is  shown  in  Figure  7.  4.  2  and  shows  that  for  a  given 

error  signal,    ^   ,   there  are  two  solutions  for  roll  rate.    The 

solution  at  the  higher  rate  is  unstable  for  the  same  reasons  that 

the  initial  point  at  p  is  unstable  in  Figure  7.  2.  4.    Therefore, 

max  *= 

to  avoid  all  possibilities  of  the  position  control  mode  being  un- 
stable the  error  signal  from  the  sensor  is  limited  to  that  giving 
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a  maximum  gimbal  angle  of  30     which  allows  the  position  control 

loop  to  drive  the  rate  of  the  spacecraft  up  to  one-half  p  .    This 

^  t^  ^  ^max 

limiting  also  enables  the  stability  analysis  of  the  position  control 
equations  to  be  accomplished  with  a  constant  gimbal  angle.    Thus 
for  small  gimbal  angles  the  position  control  equation  7.  4.  4  can 
be  written  as  follows. 


0.  192co 


</>  = 


o 


0    + 
r 


S+3 


(S+  0.  577  ojY 


(0.577.  J 


M    /I 

X       X 


(Eq.  7.  4.  10) 


For  zero  sensor  input  the  response  to  an  impulsive  moment  dis- 
turbance is  given  by 
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<^  = 


^(O. 


577go     t    +t 
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-0.  577  00    t 
e  o 


(Eq.  7.  4.  11) 
The  response  to  a  step  moment  disturbance  is  given  by 


M. 


0=' 


3  Uj  (t) 
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(0.577.  J 


-> 


t^^ 


3t 


(b.  577  00^)         (b.  577co^)' 


,577a)    t 
o   V 


(Eq.  7.  4.  12) 

For  zero  moment  disturbance  the  response  to  a  step  position 
input  angle  is  given  by 


^  =  (l>^   <Uj(t) 


(b.  577  a.  Y    „ 

-^^ ^ 2Z  t^  +0. 


577aj    t+  1 
o 


-0.  577  aj    t 
o 


(Eq.  7.  4.  13) 
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The  response  to  an  impulsive  input  is  given  by 

(Eq.  7.  4.  14) 


'0.  577u    f  <?      -0.  577u    t' 


*=*o 


The  response  is  seen  to  be  well  behaved  at  small  gimbal 
angles;  however,    consider  the  location  of  the  poles  of  the  closed 
loop  position  control  mode  as  the  gimbal  angle  becomes  larger. 
Figure  7.  4.  3  shows  the  migration  of  the  poles  for  the  value  of 
K     chosen  in  equation  7.  4.  3.    This  figure  indicates  that  at  gimbal 
angles  greater  than  30  degrees  the  system  becomes  very  oscilla- 
tory, and  at  approximately  84  degrees  the  system  becomes  unstable 
The  effect  of  the  initial  choice  in  open  loop  damping  ratio,   ^    , 
is  shown  in  Figure  7.  4.  4  which  is  plotted  for  the  same  K   .    A 

smaller  K     would  of  course  raise  the  curve  of  Figure  7.  4.  4  so 
X  ^ 

that  in  the  end  a  singular  point  at  90  degrees  would  result. 
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Figure  7.4,2      Steady  slate  roll  i-ate  for  constant 
sensor  error  signal,       ^g__    of  0,3  is  an  example  only. 
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P'igure  7.4.4    Plot  of  Boundary  between  Stable  and 
Unstable  Operation  in  Position  Control  Mode . 
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7.  5         Adaptive  Mode 

It  has  been  stated  earlier  that  the  four  gyro  system  has 
the  adaptive  feature  that  permits  continuous  three  axis  control  in 
the  event  that  one  of  the  four  gyros    is     lost.    This  mode  can  be 
illustrated  by  assuming  the  loss  of  the  number  four  gyro  control- 
ler and  examining  the  resulting  spacecraft  response.    At  the  out- 
set,   let  it  be  assumed  that  the  controller  fails  such  that  it  loses 
its  angular  momentum  at  a  steady  rate  over  time,   t..    Assuming 
no  interaxial  coupling  this  gives  a  moment  disturbance  about  the 
roll  axis  as  follows. 

•—  =   M  (Eq.    7.  5.  1) 

dt  X 

The  nnoment  disturbance  is  seen  to  be  a  step  of  H/t^  which  lasts 
for  t.   seconds,   and  this  gives  a  total  torque  impulse  of  H. 

The  spacecraft  response  to  such  a  disturbance  can  be  ob- 
tained from  the  equations  of  the  previous  section;  however  to 
illustrate  the  response  for  the  adaptive  mode  the  problem  was 
set  up  on  an  analog  computer  and  the  results  are  shown  in  Figure 
7.  5.  1  and  7.  5.  2  for  the  case  where  a  sensor  is  providing  attitude 
errors.    The  zero  input  mode  also  operates  to  provide  rate 
stabilization  upon  the  loss  of  any  controller. 
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Figure  7.5,1    Spacecraft  Roll  Rate  and  Roll  Angle  with  Gimbal  Rate 

and  Angle  for  the  Adaptive  Mode,     Data  is  from  Analog  Computer, 

I    -10^    Ib-ft-sec^,     J  =  10    Ib-ft-sec^,     t,  =  1000  sec.     H=10^ 
X  1 

Ib-ft-sec,     /,=^0.866.     See  Figure  7.5.2  for  Phase  Trajectory, 
Complete  Run  includes  Initial  Transient  and  Recovery  Transient. 
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Figure  7.5.2    Phase  Trajectory  for  Adaptive  Mode 
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7.  6  Summary 

This  chapter  has  evaluated  the  response  of  the  system 
to  input  disturbances  and  input  commands  for  the  zero  input  mode, 
the  rate  control  mode,   and  the  position  control  mode.     The  re- 
sults are  shown  in  the  illustration  of  the  location  of  the  poles  of 
the  various  modes  shown  in  Figure  7,  6,  1.     This  plot  shows  that 
an  initial  choice  of  the  damping  ratio  of  the  zero  input  mode,    £,   , 
of  approximately  0.  866  gives  satisfactory  performance  for  the 
rate  control  mode  and  the  position  control  mode;  therefore,   no 
adaptive  type  change  in  the  damping  of  the  gimbal  or  tandem,  com- 
pensation is  required  to  satisfy  all  three  modes  of  operation. 
The  effects  of  increase  in  gimbal  angle  is  seen  to  make  the  zero 
input  mode  and  the  rate  control  mode  less  oscillatory,   whereas 
the  position  control  mode  becomes  more  oscillatory.     For  a  fixed 
value  of  gain,  K    ,   there  exists  a  gimbal  angle  for  the  position 
control  mode  which  makes  the  system  unstable. 

The  system  operating  in  the  position  control  mode  has 
inherent  adaptive  characteristics  in  that  the  failure  of  any  one  of 
the  gyro  controllers  will  automatically  be  compensated  for  by  the 
repositioning  of  the  opposite  pair  of  controllers. 

The  examples  have  been  given  for  the  roll  equations. 
The  pitch  equations  are  almost  identical  except  for  the  sign  of 
the  gain  K   .     The  general  solution  for  the  yaw  equations  are 
identical  except  that  the  functions  F-      and  F^     have  been  defined 
for  a  pair  of  controllers  whereas  G     is  defined  for  a  single  con- 
tr oiler.     Also,   the  yaw  controller  uses  all  four  controllers,   and 
therefore  faster  response  is  available,   if  desired. 
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Figure  7.6,1    Illustration  of  the  location  of  the  poles  of  the  Zero 
Input  Mode         (-;„  =0.866),   the  Rate  Control  Mode  (m=l),   and 
the  Position  Control  Mode.     I'he  poles  are  shown  for  gimbal  angles 
of  zero  and  30  degrees  to  show  effects  of  change  in  gimbal  angle. 
The  open  loop  poles  of  the  position  control  mode  are  the  closed 
loop  poles  of  the  zero  input  mode. 
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CHAPTER  8 
CONCLUSIONS  AND  RECOMMENDATIONS 

8.  1 Summary  of  the  General  Concepts  of  the  Thesis 

The  thesis  is  concerned  with  the  problem  of  attitude 
control  of  a  spacecraft  engaged  in  an  extended  mission.     The 
foremost  application  of  such  a  spacecraft  is  a  manned  explora- 
tion of  the  planet  Mars.     The  factors  considered  in  choosing  a 
suitable  attitude  control  system  are 

Maximum  reliability. 

Minimum  ejection  of  mass, 

Minimum  average  power. 

Minimum  system  weight,   and 

Minimum  peak  power. 

Attitude  control  systems  which  do  not  expend  fuel  mass 
but  derive  control  moments  by  a  time  rate  of  change  in  angular 
momentum  of  a  mass  that  remains  within  the  spacecraft  are 
called  momentum  exchange  type  attitude  control  systems  and  con- 
sume only  power  in  their  operation.     In  a  space  environment, 
power  is  more  available  and  may  actually  be  re- supplied  from  the 
sun,   whereas  expulsion  control  fuel  is  limited  to  that  initially 
loaded  aboard  the  spacecraft.     Based  on  this,   the  main  theme  of 
the  thesis  is  the  study  of  momentum  exchange  type  attitude  con- 
trol systems. 

The  thesis  proceeds  to  accomplish  five  objectives  set 
forth  in  section  1.2.     In  short,   these  objectives  are 

derive  the  equations  of  motion, 

determine  the  torque  disturbances. 
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select  a  specific  control  system, 

determine  the  response  of  the  spacecraft,   and 
present  the  conclusions. 

The  attitude  control  system  chosen  to  best  satisfy  the 
spacecraft  requirements  is  a  configuration  of  four  controllers 
which  operate  to  apply  gyroscopic  torques  to  the  spacecraft. 
Each  controller  may  be  described  as  a  gyroscopic  type  device 
having  two  degrees- of- freedom.     The  four  controllers  are  ar- 
ranged in  two  pairs  with  each  pair  operating  back-to-back.     One 
pair  of  the  controllers  is  actuated  symmetrically  to  generate 
torques  applied  to  roll  the  vehicle,   whereas,   the  other  pair, 
mounted  normal  to  the  first  pair,   provides  torques  for  pitch  con- 
trol.    All  four  of  the  controllers  provide  torques  for  yaw  control. 
This  configuration  of  four  controllers  has  a  redundancy  in  its 
ability  to  effect  a  change  in  angular  momentumi  along  each  of  the 
three  spacecraft  control  axes  such  that  the  system  provides  un- 
interrupted control  upon  the  complete  loss  of  angular  momentum 
of  any  one  controller. 

The  four  controller  attitude  control  system  is  operated 
in  three  modes: 

Zero  Input  Mode 

Rate  Control  Mode 

Position  Control  Mode. 

The  zero  input  mode  enables  the  spacecraft  to  be  rate 
stabilized  in  the  absence  of  input  commands.     The  rate  control 
mode  provides  a  means  of  changing  the  attitude  of  the  vehicle  at 
the  maximum  rate  capability  of  the  controllers.     The  position 
control  mode  is  provided  for  alignment  of  the  spacecraft  with  res- 
pect to  a  reference  line  of  sight.     Satisfactory  stability  charac- 
teristics were  obtained  for  each  mode. 
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8.  2 


Summary  of  Chapter  Conclusions 


of  which 
Chapter 


Each  chapter  contains  a  brief  summary,   the  highlights 
follow. 


1:        This  chapter  introduces  the  problem  of  attitude  con- 
trol.    In  the  practical  problem  the  control  system 
engineer  must  have  knowledge  of  the  physical  charac- 
teristics of  the  spacecraft,   the  controllability  re- 
quirements,   and  the  stability  requirements. 

Chapter  2:        Chapter  2  contains  the  equations  of  motion  of  the 
spacecraft.     The  torques  generated  by  the  control 
system  depend  on  the  attitude  of  the  controller 
rotor  with  respect  to  the  spacecraft  and  the  rates 
and  accelerations  of  the  controller  rotor  with  res- 
pect to  inertial  space.     The  equations  are  presented 
in  a  manner  which  facilitates  the  evaluation  of  each 
of  these  variables. 

Chapter  3:        One  or  more  controllers  are  combined  to  form 

various  attitude  control  systems.     Control  logic  is 
required  to  give  non- interacting  control  for  the 
roll,  pitch,   and  yaw  vehicle  attitude  variables. 
A  system  of  compensation  is  devised  to  minimize 
cross- coupling. 

Chapter  4:        This  chapter  seeks  to  determine  the  torque  dis- 
turbances acting  on  a  spacecraft.     It  is  determined 
that  the  external  torque  disturbances  can  be  mini- 
mized by  design  and  operation  of  the  spacecraft. 
The  torques  which  originate  within  the  spacecraft 
from  masses  which  do  not  leave  the  system  have 
a  zero  mean  value  and  their  effects  can  be  com- 
pensated for  by  using  a  momentumi  exchange  type 
attitude  control  system. 
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Chapter  5:        This  chapter  compared  the  various  systems  and 
determined  that  the  four  gyro  controller  best  ful- 
filled the  control  system  requirements.     The  six 
gyro  control  system  is  the  second  choice.     Inertia 
reaction  wheel  systems  have  very  low  power  effi- 
ciency,  and  they  lack  the  inherent  stabilizing  char- 
acteristics of  the  gyro  controllers.     Mass  expulsion 
systems  must  be  provided  to  some  extent  to  de- 
saturate  the  momentum  exchange  system,   but  if 
used  continuously,   large  amounts  of  fuel  will  be  re- 
quired unless  a  limit  cycle  of  the  order  of  15  min- 
utes is  provided. 

Chapter  6:        The  control  loop  was  closed  using  the  four  gyro 

controller.     A  general  solution  of  the  equations  re- 
quires machine  computation  because  of  the  non- 
linearity  and  complexity  of  the  equations.     Three 
modes  of  operation  are  considered.     Zero  Input 
Mode,   Rate  Control  Mode,   and  Position  Control 
Mode. 

Chapter  7:        The  response  of  the  spacecraft  is  found  to  be  satis- 
factory for  all  modes  of  operation  using  a  fixed 
value  of  gimbal  damping.     The  zero  input  mode  and 
the  rate  control  mode  become  less  oscillatory  at 
larger  gimbal  angles,   whereas  the  position  control 
mode  becomes  more  oscillatory.     The  error  signal 
from  the  sensor  must  be  limited  to  avoid  a  possible 
unstable  condition  which  drives  the  gimbal  to  their 
full  90  degree  position. 
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8.  3 Summary  of  Characteristics  of  Gyro  Type  Controllers 

The  following  conclusions  are  presented  concerniaf^ 
gyro  type  controllers  used  to  provide  attitude  control  of  space- 
craft. 

a.  The  use  of  gyro  controllers  operating  back-to-back 
eliminates  large  cross  control  moments  typical  of  single  gyro 
controllers,   and  this  becomes  more  important  as  the  gimbal 
angle  is  increased. 

b.  Gyro  controllers  should  be  operated  without  an  ap- 
preciable amount  of  saturation  in  order  to  minimize  gyroscopic 
cross  coupling  moments.     This  statement  is  true  about  any  mo- 
mentum exchange  system,   and  suggests  that  zero  angular  mo- 
mentum type  systems  are  preferred  to  those  which  do  not  have 
zero  angular  momentum  in  their  initial  configuration. 

c.  Except  for  the  unstable,   full  90  degree  position  of  the 
gyro  controllers  where  no  control  moment  is  required,   gyro 
controllers  require  continuous  control  moments  to  provide  an 
attitude  rate  to  the  spacecraft. 

d.  During  the  angular  acceleration  of  the  spacecraft  to 
an  attitude  rate  the  torque  multiplication  for  the  gyro  torquer  is 
high.     The  ideal  torque  gain  for  a  pair  of  controllers  operating 
open  loop  is  equal  to 


1 


/= 


e.  Gyro  controllers  provide  inherent  rate  stabilization, 

f.  Since  the  gyro  gimbal  is  free  to  change  its  attitude 
with  respect  to  the  spacecraft,   generally  a  control  system  using 
several  gyro  controllers  will  have  inherent  adaptive  characteris- 
tics,  in  that,   the  failure  of  a  single  gyro  unit  will  be  compensated 
for  by  the  repositioning  of  one  or  more  of  the  other  gyros. 
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g.       In  the  absence  of  external  torques  the  spacecraft  as- 
sumes a  particular  attitude  rate  corresponding  to  a  controlled 
gimbal  position  as  given  by  the  following  equation  for  roll  rate, 

p  =  p  sin  7. 

h.       Gyro  controllers  are  more  efficient  than  wheel  con- 
trollers since  they  can  effect  a  change  in  angular  momentum  of 
the  spacecraft  without  a  change  in  their  kinetic  energy,   whereas 
wheel  controllers  require  a  kinetic  energy  change  to  effect  an 
angular  momentum  change  which  results  in  a  wheel  efficiency 
ratio  of  the  order  of  the  ratio  of  the  wheel  moment  of  inertia  to 
the  spacecraft  moment  of  inertia. 

i.       The  provision  of  a  position  control  loop  for  the  gimbal 
angle,    as  is  done  in  the  rate  control  mode,    causes  a  deterioration 
in  the  ability  of  the  gyro  to  provide  inherent  rate  stabilization. 

j.        The  ratio  of  moment  of  inertia  to  damping  coefficient 
for  the  gyro  controller  should  be  greater  than  that  typically  used 
for  integrating  gyros  of  inertial  navigation  systems 

k.       The  power  requirement  for  a  gyro  controller  is  ap- 
proximately constant  as  compared  to  an  inertia  wheel  which  has 
severe  peak  power  requirements, 

1.        Gyro  controllers  can  provide  fast,    accurate,    and  well 
damped  control  for  manned  spacecraft. 
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8,  4  Recommendations  for  Further  Study 

As  an  extension  to  the  thesis  it  is  considered  desirable 
to  continue  the  study  in  the  following  areas. 

a.  Optimization  of  the  design  of  gyro  controllers. 

b.  Study  of  controllers  which  employ  fluids  to  obtain  a 
change  in  angular  momentum. 

c.  By  the  use  of  machine  computation,    determine  the 
effects  of  the  cross  coupling  moments  acting  on  the  spacecraft. 

d.  Determine  the  feasibility  of  using  existing  angular 
moment  such  as  power  turbines  in  a  spacecraft  to  achieve  atti- 
tude control. 

e.  Investigate  devices  which  can  store  energy  as  well  as 
deliver  energy  for  use  with  inertia  reaction  wheel  control  sys- 
tems , 
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APPENDIX  A 
SYMBOLS  AND  MATRIX  NOTATION 

A.l  Matrix  Notation 

The  treatment  of  an  analytical  problem  that  uses  several 
sets  of  orthogonal  Cartesian  coordinates  is  simplified  if  matrix 
methods  are  employed  so  that  the  coordinate  frames  are  unam- 
bigously  inferred  by  the  matrix  expressions.    In  reference  1  there 
is  contained  a  particularly  good  explanation  of  matrix  notation  as 
applied  to  control  system  problems.    Thus  for  detailed  explanation 
the  reader  is  referred  to  that  reference  and  only  a  brief  discussion 
of  the  notation  is  presented  herein. 

A. 2  Coordinate  Transformation  Matrices 

The  symbol  Q  is  defined  as  an  orthogonal  transformation 
between  two  Cartesian  coordinate  frames,    and  it  is  a  3  x   3  square 
matrix.    Each  Q  shall  contain  a  double  subscript  indicating  the 
pertinent  coordinate  reference  frames  in  a  "to-from"  sequence 
from  left  to  right.    For  example,    Q,,    is  the  coordinate  transfor- 
mation which  is  post  multiplied  by  a  vector  in  the  A  frame  to 
transfornn  the  vector  to  the  I  frame.    Qat  then  is  the  coordinate 
transformation  that  transforms  a  vector  from  the  I  frame  to  the 
A  frame,   and  since  all  Q  transformations  are  orthogonal  trans- 
formations the  operation  of  inverting  a  Q  matrix  is  sinnply  that  of 
taking  the  transpose  of  the  matrix. 


Qai 


QiA 


'^=     [Qj^   ^  (Eq.    A.  2.  1) 


Choice  of  the  "to-from"  sequence  in  subs.cripts  facilitates 
multiple  transformations  so  that  the  following  is  true. 

QlR    =    "^iv  %A   ^AR  <El-   A-  2-  2) 
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It  is  sometimes  necessary  to  represent  a  coordinate  frame 
by  more  than  one  letter,    as  for  example  the  gimbal  frame  is  re- 
ferred to  as  GIM.    When  there  is  chance  of  ambiguity  in  a  subscript 
a  connma  is  used  to  separate  the  two  frames.    Thus  the  coordinate 
transformation  Q.    riM  represents  the  coordinate  transformation 
to  the  A  frame  from  the  GIM  frame.    The  coordinate  frames  used 
in  this  thesis  are  contained  in  Appendix  C. 

A.  3         Relative  Velocities 


A  velocity  is  a  vector  quantity  and  is  expressed  in  matrix 
form  by  a  column  vector.    Thus  W.|-|-,  means  the  3X1  matrix  rep- 
resenting the  angular  velocity  of  frame  E  relative  to  frame  I,    ex- 
pressed in  frame  E. 

In  matrix  equations  it  is  necessary  to  have  a  means  of 
expressing  the  operation  of  a  vector  cross  product.    Thus  consider 
a  vector 


Wjj^    =    li    +    mj    + 

nk 

(Eq. 

A. 

3.1) 

and  a  vector 

Hj^  =  aT  +  bJ  4 

Ck 

(Eq. 

A. 

32) 

i 

3 

k 

The  cross  product 

W       X   h"      = 
IR           R 

1 

m 

n 

(Refer  to  ref  29 

A 

B 

C 

p  190) 

' — 

—^ 

(Eq. 

A. 

33) 

^IR  ^   ^R    "    (mC-  nB)  i    h  (nA-lC)j  +  (IB-mA)k 

(Eq.  A.  34) 

In  matrix  form  the  result  of  the  cross  product  operation  can  be 
represented  as  a  column  matrix  such  that, 

m  C  -  nB 

IR     ^%| 


W.r^'H, 


R 


nA    -    IC 
1  B    -  mA 


(Eq.  A,  3.5) 
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W-rp  and  Ht3  respectively  gives: 


Then  upon  separation  of  the  variables  associated  with 

R 

'^  -n  nn 

0  -1 

1  0 


WiRli^Hj^ 


R 


0 


n 
■m 


A 
B 
C 


(Eq.    A.  3.6) 


Therefore  we  define  from  any  column  matrix  such  as 

1 
W. 


IR 


the  3x3  antisymmetric  matrix 


m 
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Wtp-^- 


IR 


0 

n 

■m 


~ 

n 

m 

0 

-1 

1 

0 

__ 

(Eq.    A.  3.7) 


(Eq.    A.  3.8) 


Reference  1  derives  this  definition  in  a  more  rigorous 
manner  from  the  fundamental  operation  of  differentiating  a  matrix 
equation.    Suppose 


w      =0       W      +  W 
^IR       ^RA       lA         ^AR 


(Eq.    A.  3.9) 


Differentiating  with  respect  to  time 


^IR  =  QrA  Wia  +  QrA  W:a  +  War        (Eq.   A.  3.10) 


Since   QpA   Qar  =   !>   a  matrix  Qp  a   can  be  factored  from 


equation  A.  310. 


Wffi  =  Qra 


WlA  ^  QaR  QrA  WiA 


+  W^jj       (Eq.  A.  an) 
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and  by  differentiating  each  of  the  elements  of  Q-^  ^^  and  premulti- 
plying  the  result  by  Qatd  we  find  that 

Q         Qra   is  a  3  X  3  antisymmetric  matrix 
of  the  form 

^AR  ^RA  ~  ^RA"^^^  defined  in  equation  A.  3.8. 
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A.  4         Numbering  of  Controller  Elements 

When  control  systems  are  considered  that  contain  many 
elements  some  system  must  be  devised  to  keep  track  of  the 
various  elements.    For  the  lack  of  a  better  arrangement  the  con- 
trollers have  simply  been  numbered  as  the  problems  have  been 
solved.    Therefore,   the  following  is  a  list  of  the  various  numbers 
assigned  to  the  controllers  (Gyro). 

1  Gyro  with  Spin  Reference  Axis  along  y  axis 

2  Gyro  with  Spin  Reference  Axis  along  -y  axis 

3  Gyro  with  Spin  Reference  Axis  along  x  axis 

4  Gyro  with  Spin  Reference  Axis  along  -x  axis 

5  Gyro  with  Spin  Reference  Axis  along  z  axis 

6  Gyro  with  Spin  Reference  Axis  along  -z  axis 

7  Gyro  with  Spin  Reference  Axis  in  x-y  plane 
and  rotated  120     from  Gyro  1  about   z  axis. 

8  Gyro  with  Spin  Reference  Axis  in  x-y  plane 
and  rotated  -120     from  Gyro  1  about   z  axis. 

9  Gyro  with  Spin  Reference  Axis  in  y-z   plane  and 
rotated  -a     degrees  from  -y  axis  about  x  axis 

10        Gyro  with  Spin  Reference  Axis  in  y-z   plane  and 
rotated   i-o'     degrees  from  -y  axis  about  x  axis 

A  particular  control  system  may  require  single-degree- 
of-freedom  controllers,    two-degree-of- freedom  controllers,    or 
possibly  three-degree-of- freedom  controllers  in  which  the  angular 

speed  of  the  rotor  may  be  varied  in  addition  to  the  two  degrees- 

(33) 
of- freedom  in  precession  of  the  spin  axis  of  the  controller 

Therefore  the  controllers  used  for  a  particular  system  require 
further  specification  to  fully  define  their  configuration.   Accord- 
ingly in  Appendix  G  which  presents  specific  control  systems  the- 
angles  required  to  align  the  case  of  the  controller  shown  in 
Figure  B.  9  of  Appendix  B  are  given  for  each  system.    See  also 
Appendix  C.  3  for  coordinate  transformations  for  particular  con- 
trollers. 
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A.  5         List  of  Symbols 

a  Angle  used  to  Define  Gimbal  Position  in  Case.    See  Figure 

B.  10. 

y  Angle  used  to  Define  Case  Position  Relative  to  Vehicle. 

See  Figure  B.  9.    Also  used  as  a  General  Control  System 
Variable. 

|3  Angle  used  to  Define  Rotor  Position  Relative  to  Gimbal. 

See  Figure  B.  11. 

■q   ,  r)   ,  &  r/  Represent  Primary  Control  Variables.    Also  used 

X      y  z 

as  angle  without  subscript  in  Figure  B.  3. 

'V'^  Represents  Vernal  Equinox. 

fi,aj,  &i  Angular  Rotations  of  Heliocentric  Orbital  Plane 

Reference  Frame  relative  to  Heliocentric  Inertial 
Reference  Frame. 

(j),  9  and  i//  Angular  Rotations  of  vehicle-centered  Vehicle 

Reference  frame  with  respect  to  an  Inertial 
Reference  Frame. 

X      ^      &L  fu^-f-.   Refers  to  Angular  Rotations  of  Geocentric  Orbital 
Position  Reference  Frame  relative  to  the  Geo- 
centric Orbital  Plane  Reference  Frame.    See 
Figure  B. 13. 

LAT      Latitude.    See  Figure  B.  14. 

LON      Longitude.    See  Figure  B.  14 

IJ,     and  A  Angular  Rotations  used  in  Defining  the  Geocentric 

Solar  Reference  Frame.    See  Figure  B.  15. 

A  Represents  an  Angle  when  provided  with  a  subscript.    When 

used  as  a  subscript  the  symbol  refers  to  the  Vehicle 
Principal  Axis  Reference  Frame.    See  Figure  B.  8. 

E  Energy 
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H  Angular  Momentum 

M  Moment 

P  Power 

Q  Orthogonal  Coordinate  Transformation  Matrix.    (See  Ap- 

pendix C) 

R  Resistance,    ohms 

p,  q,  and  r  Vehicle  Attitude  Rates  in  Roll,    Pitch,  and  Yaw 

W-pA         Angular  Velocity  of  Coordinate  Frame  A  relative  to 
Coordinate  Frame  I  expressed  as  a  Column  Matrix 

I  Moment  of  Inertia  of  Spacecraft 

J  Moment  of  Inertia  of  Components  of  Control  System 

n  Refers  to  number  of  columns  or  rows  of  a  matrix 

•  A  dot  over  a  symbol  indicates  a  differentiation  with  re- 

spect to  time 

■^^  Represents  the  forming  of  a  Matrix  from  a  Column  Vec- 

tor.   See  Section  A.  3. 

S  Shorthand  notation  for  the  trigonometric  sine. 

C  Shorthand  notation  for  the  trigonometric  cosine. 

Subscripts 

T  Total 

CS  Control  System 

ext  External 

VF  Refers  to  Masses  in  the  Vehicle  that  are  Rigidly  Fixed 

to  the  Spacecraft 

VM         Refers  to  Masses  in  the  Vehicle  that  Move  Relative  to  the 
Vehicle  Exclusive  of  the  Rotating  Member  of  the  Control 
System  Controller. 
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N  Represents  a  summation  symbol  for  N  Controllers 

a  Armature 

i  Summation  Index.    Also  used  in  defining  the  Heliocentric 

Orbital  Plane  Reference  Frame.    See  Figure  B.  2. 

c  Refers  to  the  Case  of  the  Controller 

g  Refers  to  the  Gimbal  of  the  Controller 

r  Refers  to  the  Rotor  of  the  Controller 

X,  y,  and  z  Generally  refers  to  the  x,  y,  and  z  direction  of  an 

Orthogonal  Coordinate  Reference  Frame. 

w  Wheel 

Z^,  YN,  and  XV       Refers  to  Angular  Rotations  of  Vehicle-Cen- 
tered Vehicle  Reference  Frame  relative  to 
Vehicle-Centered  Inertial  Reference  Frame. 
See  Figure  B.  7. 

ZV,  1;  YYjandXA      Refers  to  Angular  Rotations  of  Vehicle-Cen- 
tered Principal  Axis  Frame  relative  to 
Vehicle-Centered  Vehicle  Reference  Frame. 
See  Figure  B. 8. 

YU  and  X,  GU  Refers  to  Angular  Rotations  of  Vehicle-Cen- 

tered Gyro  Case  Axis  Reference  Frame  rela- 
tive to  Vehicle-Centered  Vehicle  Reference 
Frame.    See  Figure  B.  9. 

Zl  Refers  to  angular  rotation  of  Earth.    See  Figure  B.  12. 

See  section  B,  1  for  Symbols  used  for  Coordinate  Refer- 
ence Frames. 
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A.  5         Glossary  of  Terms 


Active  Control 


Active  control  of  a  spacecraft  is  defined  as  the  operation 
of  controlling  the  spacecraft  with  a  torque  producing   con- 
trol system  operating  with  error  sensors  in  a  closed  loop. 
The  control  system  always  consumes  energy. 

Compensation 

A  signal  provided  to  a  control  loop  which  is  proportional 
to  a  signal  of  another  loop,    and  is  used  to  minimize  the 
effects  of  an  unwanted  coupling  between  these  two  loops. 
Three  types  of  compensation  are  used  in  this  report: 
gyroscopic  coupling,    cross  control  coupling,   and  space- 
craft inertia  cross  coupling. 

Control  Logic  Matrix 

The  control  logic  matrix  is  defined  as  an  n  x   3  matrix 
which  pre -multiplies  the  primary  control  variables  to 
define  individual  signals  to  the  n-degree-of- freedom  con- 
troller. 

Control  System  Coupling  Matrix 

The  control  system  coupling  matrix  is  a  3  x   3  matrix 
which  operates  on  the  vehicle  rate  variables  resulting 
from  the  combination  of  one  or  more  controllers.    If  the 
control  system  coupling  matrix  contains  no  diagonal  terms 
then  it  is  also  the  gyroscopic  coupling  matrix. 

Control  System  Input  Matrix 

A  3  X  n  matrix  which  results  from  any  arrangement  of 
terms  of  one  or  more  controllers  in  which  the  control 
system  input  matrix  operates  on  the  control  system  input 
variables. 
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Control  System  Input  "Variables 

Any  controller  may  have  one  or  more  degrees  of  freedom 
which  may  be  considered  the  input  variables  of  the  con- 
troller,   i.  e.    with  wheels,   the  input  variables  are  angular 
acceleration  of  the  wheels.    With  pure  gyro  systems  the 
input  variables  are  precession  rates.    A  combination  of 
two  or  more  controllers  give  many  degrees  of  freedom 
which  represent  the  control  system  input  variables. 

Controller 


A  momentum  exchange  device  which  is  capable  of  applying 
control  moments  to  the  spacecraft.    A  controller  in  this 
report  is  considered  to  be  a  rotating  rigid  body.    One  or 
more  controllers  are  combined  to  form  a  complete  con- 
trol system. 

Gyroscopic  Coupling  Matrix 

Defined  as  the  control  system  coupling  matrix  with  all 
diagonal  elements  replaced  by  zeros. 

Guidance 

The  guidance  of  a  vehicle  is  defined  as  the  operation  of 
controlling  the  thrust  vectors  acting  on  a  spacecraft,  such 
that,   a  desired  trajectory  is  followed.    Guidance  of  a 
spacecraft  is  required  during  thrusting  and  possibly  during 
re-entry  if  the  aerodynamic  lift  of  the  vehicle  can  be  con- 
trolled. 

Inertial  Guidance  Measurement  Unit 

The  inertial  guidance  measurement  unit  is  defined  as  an 
assemblage  of  instrumentation  to  determine  the  specific 
force  vectors  acting  on  the  spacecraft  suitable  for  per- 
forming the  guidance  and  navigation  function. 
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Momentum  Exchange 

A  momentum  exchange  control  system  is  defined  as  a  de- 
vice which  IS  capable  of  providing  control  torques  to  a 
spacecraft  by  a  time  rate  of  change  in  angular  momentum 
within  the  device.    Consequently,   the  total  mass  and  angu- 
lar momentum  of  the  combined  vehicle  and  control  system 
remains  constant. 

Momentum  Trans fe r 

A  momentum  transfer  control  system  is  defined  as  a  de- 
vice which  is  capable  of  providing  control  torques  to.  a 
spacecraft  by  ejecting  mass  from  the  spacecraft  and 
creating  a  force- impulse  normal  to  a  lever  arm  directed 
to  the  center  of  gravity  of  the  spacecraft. 

Navigation 

Navigation  of  a  spacecraft  is  defined  as  the  operation  of 
determination  of  the  position  and  velocity  (or  position  and 
velocity  deviations),   and  computation  of  the  guidance  com- 
mands necessary  to  arrive  at  the  desired  destination. 

■Res_t  JPoint 

A  spacecraft  rest  point  is.  the  attitude  which  results  in 
zero  applied  torque,   and  at  which  attitude  the  spacecraft 
exhibits  stable  static  stability  characteristics. 

Passive  Control 

Passive  control  of  a  spacecraft  is  defined  as  the  operation 
of  controlling  the  spacecraft  purely, by  means  of  an  exist- 
ing stabilizing  torque  which  acts  on  the  spacecraft.   At- 
titude control  is  achieved  by  dissipating  energy  and  the 
system  provides  inherent  sensing. 

Pointing  Accuracy 

The  pointing  accuracy  is  defined  as  the  maximum  deviation 
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of  the  output  of  a  closed  loop  with  reference  to  a  desired 
output  when  no  disturbances  are  admitted  to  the  loop. 
Pointing  accuracy  is  a  measure  of  the  stability  of  a 
spacecraft  controlled  by  an  active  control  system. 

Primary  Control  Matrix 

The  primary  control  matrix  is  a  result  of  pre -multiplying 
the  control  logic  matrix  by  the  control  systenn  input  ma- 
trix.   The  primary  control  matrix  is  then  a  3  X  3  matrix 
and  is  diagonal  or  nearly  diagonal  for  the  range  of  control 
system  input  variables  chosen  to  define  the  control  logic 
matrix. 

Primary  Control  Variables 

Three  variables  which  give  non- interacting  (or  nearly 
non- interacting)  control  in  roll,   pitch,   and  yaw. 

Rest  Position 

For  a  spacecraft  in  an  environment  in  any  time  invariant 
situation,   there  exists  a  spacial  orientation  at  which  the 
external  moments  acting  on  the  spacecraft  exactly  cancel 
and  no  naoment  acts  on  the  vehicle. 

Semi- Passive  Control 

Semi-passive  control  of  a  spacecraft  is  defined  as  passive 
control  systems  which  increase  their  damping  properties 
by  introducing  gyroscopic  effects,    or  those  systems  which 
possess  no  static  stability  but  have  strong  damping  charac- 
teristics such  as  spin  stabilized  vehicles. 

Saturate  (and  desaturate) 

Any  momentum  exchange  type  control  system  will  have 
an  initial  angular  momentum  disposition  and  a  nnaximum 
angular  momentumi  storage  capability  for  a  particular 
control  axis.    Saturation  is  defined  as  the  percent  change 
in  the  angular  nnomentum  of  a  particular  axis  when  the 
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vehicle  is  non-rotating  with  respect  to  the  reference  frame, 
Desaturation  is  defined  as  providing  an  external  moment 
to  the  spacecraft  which  tends  to  return  the  monnentum 
exchange  control  system  to  its  initial  disposition. 

Vehicle  Attitude  Rate  Variables 

The  vehicle  attitude  rate  variables  are  the  roll,   pitch, 
and  yaw  rates  of  the  vehicle  with  respect  to  inertial 
space. 
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B.  1 

II 

B.  2 

H 

B.  3 

B 

APPENDIX  B 
COORDINATE  REFERENCE  FRAMES 

B.  1         Summary  of  Coordinate  Reference  Frames  Defined 

Figure         Symbol Title 

Heliocentric  Inertial  Reference  Frame 

Heliocentric  Orbital  Plane  Reference  Frame 

Heliocentric  Orbital  Position  Reference 
Frames 

B.  4  </>,  0,(A         Vehicle-Centered  Solar  Orbital  Reference 

Frame 

Vehicle-Centered  Inertial  Reference  Frame 

Vehicle-Centered  Velocity  Reference  Frame 

Vehicle-Centered  Vehicle  Reference  Frame 

Vehicle-Centered  Principal  Axis  Frame 

Vehicle-Centered  Gyro  Case  Axis  Reference 
Frame 

B.  10  GIM         Vehicle-Centered  Gyro  Gimbal  Axis  Reference 

Frame 

B.  11  R  Vehicle-Centered  Gyro  Rotor  Axis  Reference 

Frame 

B.  12  III  Geocentric  Inertial  Non-rotating  Reference 

Frame,   and 
E  Geocentric  Earth  Reference  Frame 

B.  13  K  Geocentric  Orbital  Plane  Reference  Frame, 

P  and  Geocentric  Orbital  Position  Reference 

Frame 
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B.  5 

I 

B.  6 

VR 

B.  7 

V 

B.  8 

A 

B.  9 

GU 

B.  14  G  Geocentric  Longitude- Latitude  Grid  Reference 

Frame 

B.  15  S  Geocentric  Solar  Reference  Franne 

B.  16  O  Vehicle -Centered  Planet  Orbital  Reference 

Franne 

An  attempt  has  been  made  to  keep  the  defined  reference 
frames  identical  with  those  given  by  Ogletree  in  reference  21. 
For  reference  frames  centered  at  a  planet  of  the  solar  system 
other  than  earth,    it  is  considered  that  Figures  B.  12  and  B.  16 
can  be  applied  with  a  suitable  subscript  denoting  the  planet  con- 
cerned. 
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HELIOCENTRIC  INERTIAL  REFERENCE  FRAME  IT 


VERNAL  EQUINOX 


FIGURE    B.  I     An  illustration  of  an  iner  tiol  fixed  frame  assuming  the  sun  as  a  fixed  point  in  space 
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HELIOCENTRIC  ORBITAL  PLANE  REFERENCE 
FRAME  H 

AXIS  Xh  passes  through  perihelion  of  spacecraft  orbit 


-ORBIT  OF  SPACECRAFT  AS 

PROJECTED  ON  CELESTIAL  SPHERE 


-VERNAL 
EQUINOX 


ASCENDING  NODE 
OF  ORBITAL  PATH 


0<i<TT 


PROJECTION  OF- 
ORBIT  OF  EARTH  ON 
CELESTIAL  SPHERE 


'H.n 


(CwCil-Sa)SilCI)(CwSil  +  SwCflCI)(Stt)SI) 
(-S«Cii-C<jSilCi)(-SwSfl.+C<JCi2cl)(SIC(u) 
SiSil  -SICil  CI 


FIGURE  B.  2  An  illustration  of  on  Orbital  Plane  Reference  Frame  defined  as  Three  Rotations  from  the 
HeiiocentriG  Inerliol  Reference  Frame, 


238 


HELIOCENTRIC  ORBITAL  POSITION  REFERENCE 
FRAME  B 


Zh.Z 


SPACECRAFT  POSITION- 


W, 


HB" 


'BH 


+Cos7;  +  Sim7    o 

—Sin  77  + Cos  •77    0 

0  0        I 


FIGURE  B.3    An  illustration  of  a  Heliocentric  Orbifol  Position Referenr.p  Frome  defined  by  a  Single 
Rotation  from  the  Heliocentric  Orbital  Plone  Reference  Frame. 
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VEHICLE-CENTERED  SOLAR  ORBITAL  REFERENCE 
FRAME  (THREE  CASES)  0.^,  ft  ^ 


PERIHELION 


^4>B 


-10  0 

0  1    0 

0+1    0 

0+1    0 

Q0B  = 

1    0  0 

v= 

0  0-1 

0   0-1 

0  0-1 

-10  0 

FIGURE  B.4  An  illustration  of  a  Vehicle-Centered  Orbital  Reference  Frame  as  related  to  the  Heliocentric 
Orbital  Position  Reference  Frame. 
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VEHICLE-CENTERED  VEHICLE  REFERENCE  FRAME  I 
SHOWN  RELATIVE  TO  VEHICLE-CENTERED  INERTIAL 
REFERENCE  FRAME  I 


MASS  CENTER  OF  VEHICLE 


INTERMEDIATE 
V  AXIS  Yy 

Tw 


1.  Frame  I  is  non-rotating  with  respect  to  fixed  stars,  i.e.  with  respect  to 

FRAME  H. 

2.  Directions  OF  Xj  ,Yi  ,Zi    are  defined  to  coincide  with  Xq.Yo  .Zq  at  time 
zero  for  any  particular  set  of  initial  conditions. 

3.  Order  of  rotations  to  place  frame  i  in  coincidence  with  frame  i. 

A.  Rotate  ABOUT  Axis  Z I  through  Angle  i/r. 

B.  Rotate  ABOUT  Axis  Y^  through  Angle  R 

C.  Rotate  about  Axis  X^  through  Angle  </>. 


FIGURE  B-5  An  illustration  of  relation  between  ttie  Vehicle-Centered  Inertial  Reference  Frame  and 
ttie  Vetiicle- Centered  Vehicle  Reference  Frame. 
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VEHICLE-CENTERED  VELOCITY-REFERENCE 
FRAME  VR 


YvR  (ALONG  VEHICLE 
INSTANTANEOUS 
VELOCITY  IN  X A  Yj 
PLANE)  ^   ^ 


WvR.,^ 


zo 


OVR.,^ 


+CAzo+SAzo    0 

-SAzo+CAzo    0 

0  0       I 


FIGURE  B.6    An  illustration  of  a  Vehicle -Centered  Velocity- Reference  Frame.  Angle  Azo  is 
rotation  about  Z     to  place  Framel  in  coincidence  with  Frame  VR. 
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VEHICLE-CENTERED  VEHICLE  REFERENCE  FRAME  3C 


TO  SUN 


DIRECTION  OF 

DIRECT  ORBIT     Yi  A^S 


MASS  CENTER  OF  VEHICLE 


INTERMEDIATE 

Y  AXIS  Y^ 


Z±  NORMAL  TO  ORBIT 
^  OF  SPACECRAFT 


Order  of  rotations  to  place  Frame  <^  in  coincidence  with  Frame  2 . 

1.  Rotate  about  Axis  Za  through  Angle  A^,^. 

2.  Rotate  ABOUT  Axis  Y^    through  Angle  Ay^. 

3.  Rotate  ABOUT  Axis  Xy    through  Angle  A^v  • 
'  See  Appendix  C  for  Qyj,. 


FIGURE  8.7  An  illustration  of  a  Vehicle -Centered  Vehicle  Reference  Frame  as  defined  by  three 
rotations  from  the  Vehicle-Centered  Orbitot  Reference  Frame. 
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VEHICLE-CENTERED  PRINCIPAL  AXIS  FRAME  A 


Azv.i 


INTERMEDIATE  Y  AXIS  Yy 
Ya 


Order  of  rotations  to  place  Frame  V  in  coincidence  with  Frame  A . 

1.  Rotate  about  Axis  Zv  through  Angle  A^v.i- 

2.  Rotate  about  Axis  Yy  through  Angle  Aw 

3.  Rotate  about  Axis  X^  through  Angle  A^^- 
See  Appendix  C   for  Qav 


FIGURE  B.8  An  lllusfration  of  a  Vehicle-Centered  Principal  Axis  Frame  as  defined  by  three  rotations 
from  the  Vehicle-Centered  Vehicle  Reference  Frame . 
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VEHICLE-CENTERED  GYRO  CASE  AXIS  REFERENCE 
FRAME  GU 


^  INTERMEDIATE 

YAXISYy 


GYRO  SRA  IS  ALONG  Yqu 
GYRO  IRA  IS  ALONG  Zqu 
GYRO    OA  IS  ALONG  Xgu 


Order  of  rotations  to  place  Frames  in  coincidence  with  Frame  GU. 

1.  Rotate  about  Axis  Zy    through  Angle  y. 

2.  Rotate  about  Axis  Yy    through  Angle  Ayu« 

3.  Rotate  about    xis  Xgj  through  Angle  Ax  qu* 
See  Appendix  C  for  Qqu  yi 


Wv.gu  '^- 


'GU 


Jcx 

0 

0 

0 

<JcY 

0 

0 

0 

Jcz 

FIGURE  B.9  An  illustration  of  a  Vehicle-Centered  Gyro  Cose  Axis  Reference  Frame  as  defined  by  three 
rotations  from  the  Vehicle-Centered  Vehicle  Reference  Frame. 
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VEHICLE-CENTERED  GYRO  GIMBAL  AXIS  REFERENCE 
FRAME  GIM 


^OIM 


^GUt'^GIM 


MASS  CENTER  OF  SPACECRAFT 


GYRO  GIMBAL   SA  IS  ALONG  Y, 


OIM 


GYRO  GIMBAL  I A  IS  ALONG  Zo,„. 
GYRO  GIMBAL  OA  IS  ALONG  Xo,„. 


OIM 


-OU 


To  PLACE  Frame  GU  in  coincidence  with  Frame  GIM. 
I.     Rotate  about  Axis  Xgy  through  Angle  a  . 


W, 


OU.GIM 


^g.GIM 


0        J,Y       0 


0       J 


gz 


Qqim.ou  - 


I  0  0 

0     +COS  a    +sin  a 
0     -sin  a    +cosa 


FIGURE  B.IO  An  lllustrotion  of  a  Vehicle -Centered  GyroGlmbal  Axis  Reference  Frame  as  defined  by  a 
single  rotofion  from  the  Vehicle-Centered  Gyro  Case  Axis  Reference  Frame. 
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VEHICLE-CENTERED  GYRO  ROTOR  AXIS  REFERENCE 
FRAME  R 


^OIM 


To  PLACE  Frame  GIM  in  coincidence  with  Frame  R 
I.     Rotate  about  Axis  Y^im  through  Angle  /3. 


W 


OIM.R    - 


'R,6IM 


0 

Jn< 

0 

0 

/3 

JrR  " 

0 

vJrY 

0 

0 

0 

0 

Jrz 

cos/ 

3  0 

-sin/3 

0 

1 

0 

+ 

sin/S 

(    0 

cos/Q 

FIGURE  B.ll  An  illustration  of  a  Vehicle- Centered  Gyro  Rotor  Axis  Reference  Frame  as  defined  by  a 
single  rotation  from  the  Vehicle -Centered  Gyro  Gimbal  Axis  Reference  Frame. 
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GEOCENTRIC  INERTIAL  NON-ROTATING  REFERENCE 
FRAME  HI  AND  GEOCENTRIC  EARTH  REFERENCE 
FRAME  E 


0°  LONGITUDE 
REFERENCE  MERIDIAN 


TO  A   ^ 

FirtST  POINT        Zl 
IN  ARIES 


EQUATORIAL 
PLANE 


'E^ 


CAzr 

SAn 

0 

0 

0 

-SAzi 

CAzi 

0 

Wm.E  = 

0 

: 

0 

0 

0 

1 

ZTTRAD 
86,164  SEC 

7.292115  X 10  ^IfAD. 

FIGURE  B-12  An  illustration  of  the  relations  between  a  Geocentric  Inertial  Non-Rotating  Reference 
Frame  and  a  Geocentric  Eortti  Reference  Frame  wliich  rotates  at  Earth  Rate. 
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GEOCENTRIC  ORBITAL  PLANE  REFERENCE  FRAME  K 
AND  GEOCENTRIC  ORBITAL  POSITION  REFERENCE 
FRAME  P 


POLAR  AXIS 


SPACECRAFT 
POSITION  IN  ORBIT 


PERIGEE 


earth's 

EQUATOR 


GROUND  TRACK 


H-P  ' H-K-^H-KP 


ASCENDING  NODE 
OF  ORBITAL  PATH 

ORBITAL  MOTION 


PRECESSION  OF  ORBITAL  PLANE 
ABOUT  earth's  POLAR  AXIS 


Xk=-(Wik+Wie)(T-To)  +  XkJo 

T    :  INITIAL  VALUE  OF  TIME 

XkJo  =  INITIAL  LONGITUDE  OF  ASCENDING    NODE  OF  ORBIT 


FIGURE  B.I3  An  illustration  of  the  relotions  between  0  Geocentric  Orbital  Plone  and  a  Geocentric 
Orbitol  Position  Reference  Frame 
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GEOCENTRIC  LONGITUDE-LATITUDE  GRID 
REFERENCE  FRAME  G 


LOCAL  MERIDIAN  PLANE 


TYPICAL  ORBITAL  PATH 
(NOT  IN  LOCAL  MERIDIAN 
PLANE, IN  GENERAL) 

Xg    5  Xp 


—  SPACECRAFT 
POSITION  IN  ORBIT 


EQUATORIAL 
PLANE 


Rp  IS  GEOCENTRIC  RADIUS 
VECTOR  TO  SPACECRAFT 

Y     IS  IN  DIRECTION  OF 
WxE       *  Rp 


^OE 


ClatClon      ClatSlon     Slat 

"Slon  Clon  0 

"SlatClon     "SlatOlon     Clat 


FIGURE  B.I4  An  Illustration  of  a  Geocentric  Longitude-Latitude  Grid  Reference  Frame. 


250 


GEOCENTRIC  SOLAR  REFERENCE  FRAME  S 


nr,E 


(TOWARD 
CENTER 
OF  SUN) 


EARTH  S  EQUATORIAL  PLANE 


PLANE  OF  THE  ECLIPTIC 


As=23°26'59" 


W 


ZTTRAD 


^|S"  365.25  DAYS 


's.m 


-s 


'AS 


^As 


H-s  =/isJo  +  Wis(T-To) 
/islo  =  INITIAL  VALUE  OF  ^g 


FIGURE  B.I5  An  Illustration  of  a  Geocentric  Solar  Reference  Frame. 
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VEHICLE-CENTERED  PLANET  ORBITAL  REFERENCE 
FRAME  0 


.w. 


IE 


W 


KP 


■  SPACECRAFT 
POSITION 
IN  ORBIT 


(PERIGEE 
OF  ORBIT) 


EARTH  S  EQUATOR 


ORBITAL  MOTION 


KO^  "KP 

5oHYp 
YoaZp 

A  A 

ZqBXp 


'OP' 


0  I  0 
0  0-1 
-10  0 


FIGURE  B.I6  An  Illustration  of  o  Vehicle-Centered  Orbital  Reference  Frame. 
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APPENDIX  C 
COORDINATE  TRANSFORMATIONS 

C.  1         Coordinate  Transformations  Between  Reference  Frames 
Relating  to  Interplanetary  Space  Analysis 

(Refer  to  Appendix  B  for  definitions  of  coordinate  frames) 


Q 


H,  II 


(CooCr^  -  SojSr^Ci)        (CojSn  +  SoaC^Ci)       SooSi 


(-SojC^- CooSl^Ci)      (-SooSn -t- CooCf^Ci)     SiCc 


SiSn 


■SiCfi 


Ci 


(Eq.    C.  1.  1) 


Q 


BH       = 


Cn 

Sn 

0 

-Sn 

Cn 

0 

0 

0 

1 

(Eq.    C.  1.  2) 


Q 


(/)B 


1 

0 

0 

0 

1 

0 

0 

0 

-1 

See  Figure  B.  4  for 
%B  ^""^  ^^B 


(Eq.    C.  1.  3) 
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(Eq.    C.  1.  9) 
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TRANSFORMATIONS  FOR  SPECIAL  CASE 
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(Eq.  C.  1.  14) 


C.  2        Coordinate  Transformations  Between  Reference  Frames 


Relating  to  Planetary  Space  Analysis 
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(Eq.  C.  2.  1) 
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(Eq.  C.  2.  3) 
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C.  3         Coordinate  Transformations  Relating  to  the  Definition  of 
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(Eq.  C.  3.  6) 


SINGLE- DEGREE- OF- FREEDOM  CONTROLLERS 
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(Eq.  C.  3.  7) 
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(Eq.  C.  3.  11) 
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APPENDIX  D 


D.  1 


Relative  Velocities 
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Assuming  Orbital  Elements  are  Constant 

{Eq.  D.  1.  1) 


r^V^ 
where  -q  =    — ^ —     of  the  order  of   1  degree/day. 

^  (Eq.  D.  1.2) 

r      =     radius  of  perihelion  of  spacecraft  orbit 

V     =     velocity  of  spacecraft  at  perihelion 

r       =     radius  of  spacecraft  from  sun 
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(Eq.  D.  1.  3) 
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(Eq.  D.  1.  4) 


(Eq.  D.  1,  5) 
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(Eq.  D.  1.  6) 


(Eq.  D.  1.  7) 


(Eq.  D.  1.  8) 


(Eq.  D.  1.  9) 
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(Eq.  D.  1.  10) 


W„r-.       =  0  Assuming  Orbital  Elements  are  Constant 
KE  ^ 


(Eq.  D.  1.  11) 
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KP 


(Eq.  D.  1.  12) 


W^p    =      0 


(Eq.  D.  1.  13) 
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(Eq.  D.  1.  14) 
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D.  2  Approximate  Relative  Velocities 

Valid  when  the  frames  nearly  coincide  or  when  angular 
velocity  is  purely  about  a  single  axis.    The  subscript  A  may  be 
substituted  for  V  for  the  case  where   OyA   =  I. 
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(Eq.    D.  2.  1) 
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which  is  further 
defined  as  a  column 
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(Eq.    D.  2.  2) 
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(Eq.    D.  2.  3) 


W 


V,  GU 


0 
0 

7 


for  the  particular  problem 
where  the  case  has  only  one 
degree  of  freedom. 


(Eq.    D.  2.  4) 
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(Eq.    D.  2.  5) 
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for  the  particular  problem  where 
j3  »  a-  and  jS  »  7    which  assumes  the 
rotor  spin  as  large. 

(Eq.    D.  2.  6) 
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D.  3  Relative  Velocities  for  Controller  1 


For  the  special  case  where  the  controller  is  in  position 


1  ,    i.  e.     A 


=  0,   A    ^^  =  0,   and  7  is  a  control  variable  , 


the  following  relative  velocities  apply, 
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(Eq.    D.  3.  1) 
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(Eq.    D.  3.  2) 
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(Eq.    D.    3.  3) 
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(Eq.    D.  3.  4) 
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(Eq.    D.  3.  5) 
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(Eq.    D,  3.  6) 
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D.  4  Time  Rate  of  Change  of  Angular  Velocities 

•       »    •    •     •  • 

For  the  case  where  f5  »  Q,y ,  p,  q  and  r.    Controller  1 

(p  (-Si3C7-SaC/3ST)+q(-Si3ST+Sc!'C/3C7)-r  CQ'Ci3-TCc^C/3-(^S/^}  3 
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W       = 


(Eq.    D.  4.  1) 


also  jS  =   0 
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APPENDIX  E 


EXACT  MOMENT  EQUATIONS  FOR 
CASE.    GIMBAL  AND  ROTOR  TERMS 


In  accordance  with  the  definition  of  the  GU,   the  GIM, 
and  the  R  coordinate  reference  frames  of  Appendix  B,   the  case 
containing  a  gimbaled    rotor  is  chosen  so  that  the  spin  reference 
axis  of  the  rotor  is  aligned  along  the  positive  y  axis  of  the  V  frame. 
This  gyro  position  is  also  called  gyro  1.    It  is  assumed  that  the 
V  frame  (vehicle- centered  vehicle  reference  frame)  and  the 
A  frame  (vehicle- centered  principle  axis  frame)  coincide. 

Equation  E.  3.  1  is  applicable  to  problems  concerning 
inertia  reaction  wheels,   but  it  is  not  particularly  suited  for 
wheels.   Accordingly,    section  F.  2  of  Appendix  F  is  a  derivation 
of  an  equation  specifically  for  inertia  reaction  wheels. 
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APPENDIX  E 
EXACT  EQUATIONS  FOR  CASE,    GIMBAL,   AND  ROTOR  TERMS 


E.  1         Case  Terms 


From  Equation  2.3.6 
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(Eq.    E.  1.  3). 
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Substituting  these  matrices  and  performing  the  indicated 
operations  gives  for  the  contribution  of  the  case  terms  the 
following  moments. 
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E.  2  Gimbal  Terms 

By  an  identical  procedure  as  the  previous  section,   the 
gimbal  terms  are  found  to  be  the  following. 
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The  preceding  is  adequate  when  two  of  the  three  attitude 
angles  are  small,    otherwise  we  must  substitute 

p  =  (/)-(//  se 

q  =  0C(/)  -  ipS<t>Cd 
r  =  i//C0C0  -  eS({) 


285 


E.  3  Rotor  Terms 

By  an  identical  procedure  as  in  section  E.  1  the  exact 
rotor  terms  for  a  rotor  of  a  controller  in  position  1  are  found 
to  be  the  following. 

Define 

J    J  -    dimetrical  moment  of  inertia  of  rotor 
rd 

J       =    polar  moment  of  inertia  of  rotor 
rp        ^ 

The  following  equations  are  adequate  when  two  of  the 
three  attitude  angles  are  small.    Otherwise  we  must  substitute 

p  =  <f)  -  ^se 

q  =  0C(^  -  JS(^C9 

r  =  (^C(^C0  -  es<^ 


286 


CM 


00 


CO 


CM  CO 


<^ 

^ 

Ul 

m 

(^ 

u 

+ 

! 

^ 

is 

U 

U 

?! 

f- 

?- 

CD 

CO 

u 

^ 

CT) 


CO 

CD 
Eh 

U 

o 


o 


a* 


^ 
?- 

?- 

T 

m 

U 

en 

•  ?- 

•^- 

1 

+ 

+ 

en 

S-. 

Sh 

a 

1 

+ 

^- 

cn 

cr 

+ 

U 
a, 


U 

u 

+ 
en 
U 


+ 

•  a 

<- 

+ 

u 

•^■ 

^ 

en 

?- 

cr 

en 

1 

ex 

c- 

1 

cn 

c- 

cj 

U 

Ul 

cr 

a 


u 

en 
cr 


U 


en 

T^ 

en 
cr 

+ 
U 


T3  a       T3  73  a        ^3 

^  Sh  ^  Sh  ^  ^H 

•-3         t-:        1-3  i-D        •-:)         i-s 


0) 
X! 


II 

si 


II 


oo 


CM  OO 


L^ 


287 


APPENDIX  F 


APPROXIMATE  EQUATIONS  FOR  CONTROLLERS 


F.  1         Arbitrary  Controller 

Using  Equation  2,  3.  8  and  assuming  the  angular  momentum 

(33) 
of  the  rotor  is  predominately  along  the  spin  axis  of  the  rotor, 

one  can  derive  the  approximate  equation 


l^csl     =^a,gimWa,gi^h  ,         ^w,.^o 

GIM 


A 


lA    "^A.  GIM^rJ 


where 


HrJ 


GIM 


Jp/3 


GIM 
(Eq.    F.  1.  1) 


(Eq.    F.  1.  2) 


Now         W 


lA 


P 

q 

r 


from  Appendix  D 


(Eq.    F.  1.  3) 


^A,  GIM      ^A,  GU  ®GU,  GIM 


(Eq.    F.  1.  4) 


and  can  be  found  from  the 
coordinate  transformations 
of  Appendix  C.  1 

^A,  GIM  "  ^GIM,  GU  ^A,  GU  ■*'  ^GU,  GIM 

(Eq.    F.  1.  5) 
and  can  likewise  be  found 
using  matrices  of  Appendix  C 
and  Appendix  D. 
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Note  that  equation  F.  1.  1  is  a  sum  of  two  parts.    The  first 
of  these  parts  will  form  the  primary  control  matrix  and  the 
second  part  will  form  the  gyroscopic  coupling  matrix  of  the  final 
equation  in  a  particular  system  using  one  or  more  of  the  control- 
lers. 

After  performing  the  indicated  operations,  the  equation  for 
a  single  controller  with  an  arbitrary  case  position  is  given  by  the 
following  equation. 
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It  is  not  considered  desirable  to  put  the  preceding  equation 
in  matrix  form  as  there  would  bo  no  simplification.    However, 
upon  use  of  the  preceding  equation  in  a  particular  controller 
position,   the  time  rates  of  the  case  angles  will  usually  vanish  and 
considerable  simplification  then  results.   At  this  point  it  is  de- 
sirable to  put  the  remaining  terms  in  matrix  form  to  assist  their 
summation  with  terms  from  other  controllers  with  the  same 
control  system  input  variables. 
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F.  2         Equations  for  Inertia  Wheels 


Unfortunately  equation  F.  1.  6  cannot  be  used  for  the  case 
of  inertia  reaction  wheels  because  it  has  been  assumed  in  that 
section  that  the  angular  velocity  of  the  controller  is  a  large  value 
and  is  constant.    Equation  E.  3.  1  can  be  used  for  the  X  wheel  and 
the  Y  wheel  or  any  other  wheel  that  lies  in  the  x-y  principle  axis 
franne  or  for  a  pure  z -aligned  wheel.    For  the  perfectly  arbitrary 
inertia  reaction  wheel  it  may  be  desirable  to  list  here  a  general 
expression  for  the  monnents.    If  we  consider  that  the  wheels  are 
rigidly  mounted  in  the  spacecraft  then  equation  2.3.8  can  be  used 
and  is  written  as  follows. 


M 


N  p 


CS  I  Lj 

-"  A       i=l 


Q 


AR 


+  WiA^QARHrJ       I 


R  R 

1 


(Eq.  F.  2.  1) 


R    i 


Or,   the  same  equation  written  with  respect  to  the  case  frame  is 

N 


M 


CS 


J 


i 


A       1=1  L 


Q 


A.GU 


r J  A,  GU'^    r  J 

GU  GU 


^WiA^QA,Gu"r|  I 


(Eq.  F.  2.  2) 


Since  the  case  holding  the  wheel  is  rigidly  attached  to  the  space- 


craft     W 


A.GU 


^=    0 


^CSI        ^Zl^A.GU     «rj        \W^1^QA,GU«rJ 

-"a       i=l  L  GU      ^  GU  I  , 


(Eq.  F.  2.  2) 
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The  xnost  difficult  part  of  this  derivation  is  that  of  finding 
H^    I  .A  considerable  simplification  of  the  algebra  results 

GU 
if  we  agree  in  this  section  to  let 


Q, 


GU,  A 


L 


e 
m 


c 

f 
n 


(Eq.  F.  2.  3) 


Where  the  elements  of  equation  F.  2.  3  are  given  by  equation  C.  1.  7 
for  the  case  where  the  vehicle  reference  frame  is  identical  with 
the  principle  axis  fraine.    Then  it  is  determined  that  the  angular 
momentum  with  respect  to  the  case  frame  is  as  follows. 


r  r  GU,  A       lA         r 

-•  GU  -■  GU  -■  GU 


%U,  R  ^GIM.  R 
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-J 


GU 


^X^ 


Jpd 


L*^xg 


j^b 


JpC 


X^ 


^x^" 

P 

■^        _ 
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• 

Jpf 

q 

+ 

Jp3 

J^n^ 

r 

0 

where 


J. 


R 


X 


0  0 


0  0 


X 


(Eq.  F.  2.  4) 


=  J. 


■J 


GU 


(Eq.  F.  2.  5) 


Note  that    H 


r  I  ='r\        QgU,  A  ^lA  ^  Jr  J        ^GU,  R  ^GIM,  R 

■^  GU  GU  GU 
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p_       — 1 
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0 

Jpd 

Jp^ 

Jpf 
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q 

-f 

Jp3 

JxS 

J^m 

Jx", 

• 

r 

(Eq 

0 
F.  2.  6) 

Following  through  the  indicated  operations  by  substituting 
the  above  matrices  into  equation  F.  2.  2  gives  the  following  equa- 
tion for  the  moments  contributed  by  the  inertia  reaction  wheels 
of  an  attitude  control  system. 
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A       i=l 
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0  f 
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e  -d 
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J  1 


P 

q 


^ia^Qa.  GU 


J^a 


Jpd 


Jxg 


J^b 


Jpe 


J^m 


•^x^ 


Jpf 


J^n 


(Eq.  F.  2.  7) 

Equation  F.  2.  7  has  been  written  in  two  parts  because  it  is 
convenient  to  lump  the  second  part  with  the  vehicle  since  that  part 
is  not  dependent  on  the  control  variable  ^.    The  second  part  relates 
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to  the  moments  contributed  to  the  vehicle  by  the  wheels  acting  as 
inert  masses  constrained  to  move  with  the  vehicle.    Therefore  the 
equation  is  written  as  the  following. 


M 


csj 
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A      i=l 
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Jp^ 
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+  Jp^ 
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whe 


re    d   =    -CA^^  ^^  St  +SA^^  ^^  ^Ayu  ^7 


(Eq.  F.  2.  8) 


^   ~      ^^X.  GU  ^^  *  ^^X.  GU  ^^YU  ^'^ 


^  ^^X.  GU  ^^YU 


from  equation  F.  2.  3  and  C.  1.  7 


For  an  X-wheel 
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Y-wheel 
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0 

Z-wheel 
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For  an  X-wheel 
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X- 
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Jp/3 
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2.  9) 
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For  a  Y-wheel 
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Y-wheel 
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1—              — 1 

0 

Jp^' 

+ 
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—                — . 

Jp/3 


Jp3 
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q 


(Eq.  F.  2.  10) 


For  a  Z-wheel 
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Z+wheel 
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""  0 

0 

+ 

Jp'^' 

-Jp(3 


JpS 


p 
q 

r 
(Eq.  F.  2.  11) 
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F.  3         Equations  for  a  Particular  Three  Degree -of- Freedom 
Controller 

Only  one  controller  is  considered  wherein  three  degrees- 

of- freedom  of  the  rotor  are  allowed.    This  controller  is  aligned 

with  its  spin  reference  axis  along  the  positive  x-axis  of  the 

vehicle- centered  principal  axis  frame.    Thus  in  equation  E.  3.  1 

if  we  let    7  -  7- 00     and  .T        =  0,    the  following  equation  is  obtained 
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-00-37 
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CaCy 
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(Eq.  F.  3.  1) 


298 


F.  4         Equations  for  Particular  Two-Degree-of- Freedom 
Controllers 

In  the  synthesis  of  the  nnoment  equations  for  those  control 
systems  that  are  composed  of  several  controllers  of  the  two- 
degree-of- freedom  type,    it  is  convenient  to  have  the  simplified 
monaent  equations  of  each  controller.    Thus  the  controllers  de- 
fined by  the  following  angles  are  presented.    See  Figure  B.  9  and 
paragraphs  A-4  and  C-3. 

TABLE  F.  4 
POSITION  ANGLES  FOR  TDF  CONTROLLERS 


TDF  CONTROLLER 

CASE  ANGLES                              | 

^X,  GU 

^YU 

7 

.  1 
2 
3 
4 
7 
8 

0 
0 
"0 
0 
0 
0 

0 
0 
0 
0 
0 
0 

0     +7 

±180°  +  7 

-  90°  +  7 

+  90°  +  7 

+120°  +  7 

-  120°  +7 

If  the  above  reference  angles  are  substituted  into  equation 
F.  1.  6  the  following  equations  are  obtained.    Note  that  7  is  a  case 
angle  retained  as  a  control  input  variable. 
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F.  5 Equations  for  Particular  Single -Degree -of- Freedom 

Controllers 

In  the  synthesis  of  the  moment  equations  for  those 
control  systems  that  are  composed  of  several  controllers  of  the 
single- degree- of- freedom  type,    it  is  convenient  to  have  avail- 
able the  simplified  equations  of  each  controller.    The  following 
table  of  angles  define  the  orientation  of  the  case  relative  to  the 
vehicle  centered  principal  axis  frame.    See  Figure  B.  9, 

TABLE  F.  5 
POSITION  ANGLES  FOR  SDF  CONTROLLERS 


SDF  CONTROLLER 

•  CASE  ANGLES 

^X,  GU 

^YU 

T 

1 

0 

0 

0 

2 

0 

0 

±180° 

3 

0 

-90° 

-90° 

4 

0 

+90° 

+90° 

5 

+90° 

0 

-90° 

6 

-90° 

0 

+90° 

7 

0 

-90° 

+  120° 

8 

0 

-90° 

-120° 

9 

-a 
o 

0 

±180° 

10 

o 

0 

±180° 
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SINGLE-DEGREE- OF- FREEDOM  CONTROLLERS 
(Subscripts  of  Variables  Omitted)-'' 


Controller  1 
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SDFll 


J£3 


A 


qSc    -  r  Cq'" 

0 

- aSa    -  pSa 

+  aCa  +  pCa_ 


(Eq.  F.  5.  1) 


Controller  2 
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(Eq.  F.  5.  2) 
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(Eq.  F.  5.  3) 
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(Eq.  F.  5.  4) 
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=  Jp0 

-p  Ca  -    r  So- 
_(q  -  a)  S  a  _ 

(Eq.  F.  5.  5) 
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(Eq.  F.  5.  G) 
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(Eq.  F.  5.  7) 
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(Eq.  F.  5.  8) 
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(Eq.  F.  5.  9) 
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(Eq.  F.  5.  10) 
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'!'    Note  that  the  equations  of  this  section  have  been  written  in 
the  form  of  a  column  vector  because  they  are  simple  expres- 
sions with  only  a  few  terms.    When  a  number  of  the  above 
equations  are  added  to  form  a  complete  control  system,   the 
form  of  Equation  3.  2.  3  is  recommended. 
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APPENDIX  G 

EQUATIONS  FOR  PARTICULAR  MOMENTUM  EXCHANGE 
ATTITUDE  CONTROL  SYSTEMS 

G.  1         INTRODUCTION 


A  spacecraft  attitude  control  system  of  the  momentum 
exchange  type  will  usually  contain  a  number  of  momentum  ele- 
ments either  of  the  inertia  reaction  wheel  or  of  the  gyro  type. 
In  this  thesis  inertia  reaction  wheels  are  referred  to  as  simply 
inertia  wheels  and  have  fixed  axes  of  rotation  whereas  gyros 
are  called  controllers  to  differentiate  from  the  conventional 
gyros  used  in  inertial  reference  systems  and  in  addition  to  its 
spin  motion  a  controller  may  have  one  or  more  degrees  of 
freedom.    To  determine  the  moment  contributions  of  a  particular 
system  consisting  of  several  controllers  or  wheels  as  the  case 
may  be,    one  needs  only  to  sum  the  contributions  of  the  individ- 
ual controllers  and  wheels.   A  number  of  the  more  common 
controllers  are  represented  in  Appendix  F. 

In  this  Appendix  a  number  of  control  systems  have  been 
defined  by  choosing  various  configurations  of  controllers.    The 
following  table  indicates  the  controllers  which  constitute  the 
system. 
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'^     The  numbers  in  the  parenthesis  partially  define  the  control 
system  by  giving  the  controllers  used  for  roll,   pitch,   and 
yaw  respectively.    Thus  the  code  following  the  Sun  Pointing 
System  (0-34-34)  indicates  that  there  is  no  roll  control; 
controllers  3  and  4  control  both  pitch  and  yaw.   A  controller 
used  in  only  one  axis  indicates  a  single  degree-of- freedom 
controller,   and  for  controllers  1  through  6  a  controller  used 
for  two  axes  indicates  a  two-degree-of- freedom  controllers. 
A  controller  used  for  three  axes  is  a  three-degree-of- freedom 
controller. 

#     These  systems  are  adaptive  versions  of  the  (12-34-1234) 
four  controller  system. 
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Single  Controller  System  (3-3-3) 


Description  of  System: 


Control  Logic: 


System  consists  of  a  number  3  con- 
troller which  has  three  degrees  of 
freedom.    The  rotor  of  the  controller 
is  free  to  precess  in  two  directions 
plus  the  rotor  is  capable  of  being  ac- 
celerated about  its  axis  of  symmetry. 

No  control  logic  is  required  for  the 
single  controller  system  since  the 
system  contains  three  degrees  of  free- 
dom, and  the  control  system  input  ma- 
trix can  be  nearly  diagonalized  by 
suitable  arrangement  of  the  matrix. 

+  j3  gives  roll  to  left 

+  7  gives  pitch  nose  down 

+  a  gives  yaw  to  left 


Moment  Contribution  of  Control  System: 


^(3-3-3) 


'JA 


L+CaST  -Cc^C7 


'+Cc^Ct 

-STCa 

-C7SQ'' 

i3 

=  jp 

^SyCa 

+Cq'C7 

-S7SCV 

/37 

L 

+s» 

0 

+  Cq'_ 

•  • 

0 

+  Sq' 

-Cq'S7~ 

[p- 

+  Jp/3 

-Sc^ 

0 

+  Cc^C7 

i 

0     J      i_rj 


(Eq.  G.  2.  1) 
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Sun  Pointing  System    (0-34-34) 

Description  of  System:       This  system  is  a  special  case  of  the 

four  controller  system.    The  controllers 
are  TDF. 


Control  Logic: 


y 


r     ~\ 

^3 

^4 

_ 

^3 

^4 

L.      _ 

0 

-1 

0 

0 

+  1 

0 

0 

0 

-1 

0 

0 

-1 

+  7o    gives  pitch  nose  down 
+  «      gives  yaw  left 


Moment  Contribution  of  Control  System: 


M 


(0-34-34) 


A 


-2Jp3 


0         -Co'CTg         +Sq'ST3 


-Cq 


t 

^x 

• 

"y 

» 

L^d 

+  2Jp^ 


-Sa 


+  30'         -Cc^St, 


+  CaST3  0 


P 

q 


(Eq.  G.  2.  2) 
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The  equation  written  in  terms  of  control  input  variables 


IS 


-] 


M 


(0-34-34)1  ^  ~  ^-^P^ 


0         +CaCy^       -SaSTg 


+  2Jp3 


0 

0  +Sa 

-Sa  0 


+  Co'ST3         ° 


+  Ca 


o1 


T. 


a 


— 

p       -, 

C»S7^ 

p 

0 

q 

0 

r 

_ 

_ 

(Eq.  G.  2.  3) 
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Sun  Pointing  System     (12-0-12) 

Description  of  Systenn:       This  system  is  a  special  case  of  the 

four  controller  system.    The  controllers 
are  TDF. 


Control  Logic: 


r-       -1 

^1 

^2 

^1 

^2 

L.     — 1 

+1 

0 

-1 
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0 

0 

0 

0 

0 

^x 

0 

^y 

-1 

\ 

-1 

+  7-1   gives  roll  right 
+  a.  gives  yaw  left 


Moment  Contribution  of  Control  System: 


M, 


'^A 


-Cq'Ct-, 

0 

-SaSy^ 

2Jp/3 

0 

0 

0 

\ 

0 

0 

-Ca 

• 

- 

L  ^ 

r 

-| 

f    1 

0 

s& 

0 

P 

+  2Jp^ 

-Sa 

0 

-CaSy^ 

q 

0 

Cc^St^ 

0 

r 

- 

L  _j 

(I 

:q.c 

5.  2 

.  A 

0 
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The  equation  written  in  terms  of  input  control  variables  is 


^(12-0-12)1  ^  ~   ^'^P^ 


0  0  0 


0 


+  2Jp3 


■Sa 


0 


Ca 


Sa 


0  CaSy, 


\ 

0 

a 

_     _ 

— , 

"    - 

0 

P 

Cq'St. 

q 

0 

r 

(Eq.  G.  2.  5) 
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Four  Controller  System  (12-34-1234) 


Description  of  System: 


Four  TDF  Controllers  are  mounted 
in  positions  1,  2,  3,    and  4. 


Control  Logic: 


-           -^ 

^1 

'^2 

^3 

^4 



^1 

^2 

^3 

^4 

■»               — 

+  1 
-1 

0 
0 
0 
0 
0 
0 


0 
0 

-1 

+  1 

0 
0 
0 
0 


'x    ! 


y 


r]      1 


+  Ti   gives  roll  to  right. 

+  Tq  gives  pitch  nose  down. 

+  a    gives  yaw  to  left. 

Moment  Contribution  of  Control  System: 

CaCTi  0 


^(12-34-1234)]^  ~  ^^P'^ 


+  2Jpi3 


0 
0 

0 
2Sq' 


0 
0 


+  C0ST3         +CaST^ 


-ScvSt 


-CaCy^        +Sq'ST3 


2Cq' 


■Cq'St. 


• 

"x 

"y 

e 

-CccSt 


(Eq.  G.  2.  10) 
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The  equation  written  in  terms  of  input  control  variables  is 


I 


^(12-34-1234)J^  -  ^"^P^ 


2Jp/3 


0  Cc^Ct, 
0  0 

0  2Sa 
-2Sa  0 


CaSTg  Cq'St^ 


Sq'STi 

"^1 

SttSTg 

^3 

2Ccv 

« 

""" 

^    - 

CO'STg 

P 

CcfST^ 

q 

0 

r 

„ 

(Eq.G.  2.  11) 
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Six  Controller  System    (56-34-12) 


Description  of  System: 
Control  Logic: 


Six  SDF  Controllers  are  mounter] 
in  positions  1,  2,  3,  4,  5,  and  6. 


_         — 

^1 

^2 

^3 

= 

^4 

^5 

^6 

_      _i 

0 

0 

-1 

0 

0 

-1 

0 

-1 

0 

0 

-1 

0 

+  1 

0 

0 

+  1 

0 

0 

y 


+  <:V(-  gives  roll  to  right 

+  Q'„  gives  pitch  nose  down 

+  0'    gives  yaw  left 

Moment  Contribution  of  Control  System: 


^(56-34-12)1^  ~  ^'^p'^ 


^^-5 

0 

0 

0 

-Cc.3 

0 

0 

0 

-Ca 

+  2Jpi3 


0 

-Sa^ 
+  Sa. 


+  Sa 

0 

+  Sar 


1 


~,   "" 

"x 

• 

"y 

« 

L 

\ 

■Sa 


■Sa 


-< 

~ 

3 

P 

5 

q 

r 

(Eq.G.  2.  12) 
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The  equation  written  in  terms  of  control  system  input  variables 
is 


'^(56-34-12)J^  =  2Jp|3 


Ca^ 

0 

0 

0 

+  C.3 

0 

0 

0 

Cg 

"■    # 

^5 

# 

^3 

9 

^1 

_       — 

+  2Jp/3 


0 
■Sa 


1 


Sq-i  -Sq-o 


0  -Sor 


3  5 


(Eq.  G.  2.  13) 
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Three  Controller  Orthogonal  System    (5-3-1) 


Description  of  System: 


Control  Logic: 


System  consists  of  SDF  Controllers 
Numbers  1,  3,  and  5. 

None  required  since  system  has 
only  three  degrees  of  freedom  and 
control  system  input  matrix  can  be 
nearly  diagonalized  by  arrangement. 

+  0-      gives  roll  to  right. 

* 
+  ao    gives  pitch  nose  down. 

+  a.     gives  yaw  to  left. 


Moment  Contribution  of  Control  System: 


^5-3-l)jA  =  Jp^ 
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q 


(Eq.  G.  2.  14) 
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Two  Controller  Yaw  System    (0-0-12) 
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(0-0-12)1^  "  ^"^P^ 


'i 


pSo' 
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c,  C  a^ 


Two  Controller  Pitch  System    (0-34-0) 
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(0-34 


00^ 


=  2J^^ 
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Two  Controller  Roll  System    (56-0-0) 


M 


(56-0-0)J^  -  ^-^P 


=  2J^/3 


-a^Ca^ 

-rsa^ 

qSa^ 

Logic 


^1  =^2 


a,   gives  yaw  to 
left 


(Eq.G.  2.  17) 


Logic 


^3  =^4 


+  0-0  gives  pitch 
nose -down 


(Eq.  G.  2.  18) 


Logic 


^5  =  ^6 


+  Q'c  =  gives  roll  to 
right 


(Eq.G.  2.  19) 
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MIT  "Vertical  Vee"    (0-0-9,  10) 


Description  of  System:  System  consists  of  two  SDF  Controllers, 


Control  Logic: 


None  proposed  since  system  is  to 
be  used  as  passive  control  system. 
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(Eq.  G.  2.  20) 
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Four  Controller  Sun  Pointing  System    (0-34-12). 


Description  of  System: 
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(0-34-12) 


J  A 


=  2Jp^ 


Six  Controller  system  without 
Controllers  Number  5  and  6. 
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(Eq.  G.  2.  21) 


Logic 


(—        -^ 

^1 

^2 

= 

• 
^3 

^4 

L      ^ 

0-10 
0-10 
0  0+1 


0 


0        +1 


y 


(Eq.  G.  2.  22) 
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Three  Orthogonal  Inertia  Wheels  (Not  a  Gyro  System) 


Control  Logic: 


None  required  because  system  has  three 
degrees  of  freedom  and  control  system 
input  matrix   can  be  diagonalized  by  inspec- 
tion. 
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